AD- A 196  221 


UNLIMITED  DISTRIBUTION 


1+ 


M  FILE  C0P> 


National  Defence  Defense  nationale 


Research  and 
Development  Branch 


Bureau  de  recherche 
et  d6veloppement 


TECHNICAL  MEMORANDUM  88/203 

February  1988 


A  SURFACE  PANEL  METHOD  FOR  THE 
CALCULATION  OF  ADDED  MASS  MATRICES 
FOR  FINITE  ELEMENT  MODELS 

T . A .  Vernon  -  B.  Bara  -  D.  Hally 


DTIC 

jELECTEI 
JUN  2  2  1988 ! 

*  i 

H 


I 


Defence 

Research 

Establishment 

Atlantic 


Centre  de 
Recherches  pour  la 
Defense 
Atlantique 


Canada 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  release; 

Distribution  Unlimited 


•'  iVf 


UNLIMITED  DISTRIBUTION 


1+ 


National  Defence 

Reaearch  and 
Development  Branch 


Defense  nationale 

Bureau  de  recherche 
et  ddveloppement 


A  SURFACE  PANEL  METHOD  FOR 
THE  CALCULATION  OF  ADDED  MASS  MATRICES 
FOR  FINITE  ELEMENT  MODELS 

T . A .  Vernon  -  B.  Bara  -  D.  Hally 


February  1988 


Approved  by  L.J.  Legget 

DISTRIBUTION  APPROVED  BY 


TECHNICAL  MEMORANDUM  88/203 


Defence 

Research 

Establishment 

Atlantic 


Centre  de 
Recherches  pour  la 
Defense 
Atlantique 


CanadS 


Abstract 


A  method  of  generating  the  added  fluid  mass  matrix  for  use  in  dynamic  analyses  of 
submerged  structures  is  presented.  This  method  uses  a  distribution  of  panel  singularities 
in  the  form  of  sources  or  dipoles  on  the  fluid/structure  interface  to  represent  the  velocity 
potential  in  the  fluid  surrounding  the  structure.  The  fluid  added  mass  matrix  is  calculated 
by  relating  the  pressure  field  in  the  fluid  to  the  structure  surface  accelerations  via  the  Euler 
equation  and  the  assumptions  of  potential  flow  theory.  Based  on  the  results  of  a  variety  of 
investigations  of  computational  performance,  convergence  and  dynamic  response  involving 
floating  cylinders,  ship  hulls,  and  propeller  blades,  the  surface  panel  method  is  shown  to 
provide  an  effective  alternative  to  the  finite  element  added  mass  matrix  system.  The  surface 
panel  program  suite  has  been  developed  in  a  form  compatible  with  an  existing  finite  element 
analysis  package,  although  the  dependence  on  the  particular  solution  system  is  not  large. 


Resume 

On  pr4sente  une  m6thode  de  generation  de  la  matrice  des  masses 
fluides  ajoutdes  utilisee  dans  1' analyse  dynamique  des  structures 
submerges.  Dans  cette  methode,  une  distribution  de  singularity  de 
plaque  sous  la  forme  de  sources  ou  de  dipSles  k  1' interface  du  fluide  et 
de  la  structure  represents  le  potentiel  des  vitesses  du  fluide  entourant 
la  structure.  La  matrice  des  masses  fluides  ajout6es  se  calcule  par 
association  du  champ  de  pression  du  fluide  aux  accelerations  k  la  surface 
de  la  structure  au  moyen  de  1' equation  d* Euler,  dans  le  cadre  des 
hypotheses  de  la  th4orie  des  ecoulements  potentials.  En  se  fondant  sur 
les  resultats  de  diverses  etudes  de  performance  de  calcul,  de  convergence 
et  de  r4ponse  dynamique  mendes  k  l'aide  de  cylindres  flottants,  de  coques 
de  navire  et  de  pales  d'heiice,  on  montre  que  la  methods  des  plaques  de 
surface  est  un  substitut  efficace  k  la  matrice  des  masses  ajoutees  k 
elements  finis.  Une  serie  de  programmes  de  plaques  de  surface  a  6t4 
eiabor4e  sous  une  forme  compatible  avec  un  actual  ensemble  de  programmes 
d' analyse  k  elements  finis,  bien  que  la  ddpendance  sur  le  systems  des 
solutions  particulibres  ne  soit  pas  grande. 
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Added  mass  in  sway 
Added  mass  in  surge 

Separation  of  field  and  source  point,  radius  of  circle 
Separation  of  imaged  field  and  source  point 
Time 

Free  stream  velocity  (rigid  body  velocity)  in  x  direction 
Characteristic  free  stream  velocity  in  x  direction 
Free  stream  velocity  (rigid  body  velocity)  in  y  direction 
Free  stream  velocity  (rigid  body  velocity)  in  z  direction 
field  point  coordinate  vector  (bold  face) 
structure  surface  point  coordinate  vector  (bold  face) 
local  x  coordinate  in  panel  coordinate  system 
local  y  coordinate  in  panel  coordinate  system 
Height  of  free  surface 


T  Structure  surface  boundary 

{5}  Nodal  displacement  vector 

fj  Surface  wave  amplitude 

p  Density 

a  Source  strength  (assumed  constant  on  a  panel) 

4>  Velocity  potential  function 

w  Characteristic  circular  frequency 


1  Introduction 


The  prediction  of  the  dynamic  structural  response  of  submerged  structures  requires  a 
knowledge  of  the  effect  of  the  fluid  on  the  structure.  In  the  general  case,  the  interaction 
leads  to  a  complex  coupled  boundary  value  problem  which  is  often  of  infinite  extent.  By 
applying  certain  assumptions  justifiable  for  the  lower  vibration  modes  of  most  structures  of 
interest,  the  interaction  can  often  be  uncoupled.  In  the  simplest  representation,  the  effect 
of  the  fluid  on  the  structure  will  be  proportional  to  only  the  interface  accelerations,  and 
thus  can  be  included  as  an  added  mass  component  in  the  dynamic  equilibrium  equations 
of  the  structure.  The  added  mass  concept  is  exactly  equivalent  to  an  integral  of  the  total 
kinetic  energy  which  the  body  motion  imparts  to  the  surrounding  fluid  in  a  potential  flow 
solution.  For  the  numerical  solution  of  the  equilibrium  equations  by  finite  element  methods, 
the  added  mass  is  in  the  form  of  a  fully  populated  mass  matrix  which  is  combined  with  the 
structural  mass  matrix. 

Several  methods  are  available  for  the  calculation  of  the  added  fluid  mass  matrix  for  finite 
element  models.  The  most  widely  used  is  probably  the  finite  element  method  itself,  whereby 
the  fluid  surrounding  the  structure  is  discretized  with  finite  elements  similar  to  those  used  in 
solid  element  modelling,  but  usually  based  on  a  single  variable  per  node,  typically  pressure 
or  velocity  potential.  This  approach  has  advantages  in  terms  of  compatibility,  and  as  a 
variationally-based  formulation,  automatically  provides  a  symmetric  mass  matrix.  The 
symmetric  form  of  both  the  mass  and  stiffness  matrices  is  important,  since  it  is  exploited 
in  the  matrix  algorithms  of  most  finite  element  solver  systems.  The  major  disadvantages 
of  finite  element  fluid  modelling  are  the  effort  required  in  the  generation  of  an  effective 
fluid  model  and  the  larger  CPU  time  and  storage  requirements  necessitated  by  these  large 
models.  As  well,  for  problems  of  infinite  extent,  the  imposition  of  the  far-field  boundary 
condition  at  a  finite  distance  in  the  finite  element  fluid  model  can  adversely  affect  the 
solution.  The  use  of  infinite  elements  which  incorporate  singular  shape  functions  provides  a 
better  representation  of  infinite  domains,  but  the  best  accuracy  is  still  obtained  when  these 
elements  are  combined  with  at  least  one  layer  of  standard  fluid  elements,  necessitating  again 
a  fluid  model. 

Since  the  interface  surface  is  the  only  part  of  the  structure  in  contact  with  a  surrounding 
fluid,  it  is  natural  to  try  to  formulate  the  interaction  problem  in  terms  of  a  boundary 
integral  approach.  This  is  possible  using  distributions  of  singularities  on  the  body  surface 
to  represent  the  potential  (or  pressure)  in  the  fluid  surrounding  the  body.  For  closed 
structures,  the  boundary  integral  approach  embodies  assumptions  about  the  fluid  potential 
on  either  side  of  the  interface  surface.  Different  assumptions  about  the  fluid  potential  inside 
the  body  lead  to  different  formulations  of  the  singularity-based  methods,  typically  cast  in 
discrete  form  as  source  or  dipole  distributions  on  panels  on  the  interface  surface.  In  the 
terminology  of  boundary  integral  methods,  both  types  of  singularity  distributions  can  be 


termed  indirect  or  direct  formulations,  depending  on  whether  the  governing  equations  are 
written  in  terms  of  velocity  potential  or  pressure  respectively.  In  many  other  applications 
of  boundary  integral  methods,  the  direct  formulation  tends  to  be  more  efficient;  however, 
both  approaches  involve  the  same  amount  of  computation  in  the  calculation  of  added  mass 
matrices.  An  example  of  the  interpretation  of  the  source  formulation  as  a  direct  method 
can  be  found  in  Reference  1. 

The  method  discussed  here  uses  the  well-known  surface  panel  singularity  distribution 
method  originally  presented  by  Hess  and  Smith2,  and  since  applied  successfully  in  the  so¬ 
lution  of  many  problems  in  the  potential  flow  regime.  The  necessary  discretization  of  the 
body  surface  into  quadrilateral  panels  lends  itself  well  to  use  with  finite  element  models 
where  the  fluid/structure  interface  elements  can  be  interpreted  as  hydrodynamic  panels. 
This  approach  has  been  discussed  by  Vorus  and  Hylarides,3’4  and  the  theoretical  develop¬ 
ment  reported  in  this  study  follows  their  work  closely.  In  addition,  a  simple  reformulation 
of  the  integral  equation  in  terms  of  a  dipole  distribution  in  place  of  the  Bource  distribution 
which  forms  the  basis  of  their  approach  is  shown  to  offer  significant  improvements  in  accu¬ 
racy  in  certain  cases,  with  no  increase  in  computing  time  demand.  As  well,  there  do  not 
appear  to  be  any  practical  evaluations  or  comparisons  of  the  results  of  the  panel  methods  to 
experimental  or  finite  element  based  mass  matrix  results  in  the  literature.  The  objective  of 
the  current  work  has  been  to  develop  the  surface  panel  method  as  a  practical  alternative  to 
an  existing  finite  element  based  fluid  modelling  capability,  and  to  evaluate  its  performance 
more  extensively  by  comparison  with  that  method. 

The  surface  panel  method  offers  several  important  advantages  over  the  finite  element 
method  in  the  generation  of  fluid  mass  matrices.  First,  fluid  modelling  is  not  required,  which 
eliminates  a  large  amount  of  data  preparation  and  input;  second,  the  formulation  is  much 
less  CPU  intensive  and  requires  less  file  storage  space;  and  third,  the  formulation  is  less 
susceptible  to  problems  of  element  aspect  ratio  and  size  disparity.  The  disadvantages  of  the 
surface  panel  method  are  largely  implementational  in  nature  rather  than  theoretical,  since 
effective  structural  modelling  does  not  always  provide  adequate  hydrodynamic  modelling. 
In  certain  cases,  problem  specific  coordinate  and  grid  transformations  are  required.  The 
surface  panel  method  also  produces  an  unsymmetric  fluid  mass  matrix,  which  can  affect  the 
overall  solution  produced  by  a  general  purpose  finite  element  solver  system  which  exploits 
an  assumed  mass  and  stiffness  matrix  symmetry.  The  evaluation  of  the  integrals  in  any 
boundary  integral  formulation  also  requires  more  effort  than  is  typically  necessary  in  a 
standard  finite  element  approach. 

Despite  the  apparent  disadvantages,  the  potential  savings  in  CPU  time  and  the  elimi¬ 
nation  of  fluid  modelling  make  the  surface  panel  method  extremely  attractive.  The  report 
presents  a  theoretical  outline  for  the  surface  panel  method,  and  discusses  its  implementation 
as  a  direct  substitute  for  the  finite  element  fluid  mass  matrix  routines  within  the  existing 
finite  element  solution  system  VAST5.  The  method  has  been  applied  in  the  analysis  of 
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several  structures  for  which  exact  or  comparison  solutions  are  available.  Rigid  body  added 
mass  predictions  are  compared  to  an  exact  solution  for  a  cylinder,  and  strip  theory  solu¬ 
tions  for  two  frigate  hulls.  Natural  frequency  predictions  using  panel  method  fluid  mass 
matrices  are  compared  to  results  utilizing  finite  element  based  fluid  mass  matrices  and  to 
experimental  results  for  a  floating  half-cylinder  and  for  two  propeller  blades.  The  computer 
processing  times  for  the  fluid  mass  matrix  generation  for  the  two  methods  are  also  reported. 
Details  of  several  of  the  major  simplifying  assumptions  and  the  integral  evaluation  methods 
used  in  the  surface  panel  formulation  are  presented  as  appendices. 

2  Theory 

A  structural  dynamic  response  analysis  by  finite  element  method  requires  the  solution 
of  the  structural  equilibrium  equations  in  the  form 

[M]«  +  [C]5  +  [*]«  =  F,(t)  +  F.(t)  (1) 

In  the  presence  of  a  surrounding  fluid,  the  effect  of  the  fluid  on  the  structure  enters 
through  the  forcing  function  F/(t),  which  ia  general  embodies  terms  proportional  to  the 
structure’s  surface  displacement,  velocity,  and  acceleration.  If  fluid  compressibility  or  a 
free  surface  is  included,  these  terms  become  frequency  dependent  and  without  simplifying 
assumptions,  the  interaction  is  of  a  complex  nature.  For  large  structures,  the  vibration 
frequencies  of  interest  are  infra-sonic,  and  the  effects  of  fluid  compressibility  are  generally 
neglected.  This  frequency  range,  on  the  other  hand,  is  often  above  the  range  in  which  free 
surface  effects  are  significant.  (A  more  detailed  discussion  of  the  assumptions  associated 
with  the  free  surface  boundary  condition  is  given  in  Appendix  A) .  If  damping  and  displace¬ 
ment  proportional  terms  can  also  be  neglected  on  the  basis  of  order  of  magnitude,  then  the 
fluid  forcing  reduces  to  a  pressure  term  proportional  only  to  the  surface  accelerations  of  the 
structure.  This  extensive  series  of  assumptions  is  commonly  made  in  the  dynamic  analysis 
of  large  submerged  structures,  and  is  usually  necessary  to  render  such  analyses  tractable. 

To  obtain  the  fluid  forces  on  the  structure,  the  relationship  between  the  fluid  pressure 
field  and  the  interface  accelerations  must  be  defined.  This  relationship  can  be  determined 
via  potential  flow  theory,  in  which  a  velocity  potential  function  ^(x)  satisfying  the  Laplace 
equation  and  appropriate  boundary  conditions  is  used  to  define  the  fluid  flow  field.  For  a 
structure  vibrating  in  an  infinite  fluid  domain,  these  boundary  conditions  are 

<j>  — ►  0  as  x  — ♦  oo  (2) 
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where  T  is  the  fluid/structure  interface  and  the  substantive  derivative  is  used  on  the  right- 
hand  side  of  (3).  If  a  free  surface  is  present 


<P<f>  d<f>  n 


(<) 


is  the  linearized  free  surface  boundary  condition  to  be  satisfied. 

A  number  of  techniques  can  be  employed  to  solve  for  the  potential  function  which 
satisfies  the  above  relations.  In  the  surface  panel  method,  function  ^>(x)  is  defined  by  a 
distribution  of  singularities  of  constant  strength  on  panels  on  the  interface  surface.  If  a 
simple  source  distribution  is  used,  the  field  potential  at  x  =  (z,  y,  z)  is  written  in  the  form 


<£(x)  =  fr  G(x,x')<r(x')  dr(x')  (5) 

where  x1  represents  the  coordinate  vector  of  a  point  on  the  panel  surface  and  a  is  the  source 
strength.  In  discrete  form,  (5)  becomes 


Np  . 

#*)  =  I  rfT(xj) 

i-i 


(6) 


where  the  constant  panel  source  strength  has  been  removed  from  the  panel  integral.  The 
kernel  function  G  in  the  current  formulation  is  the  free  space  Green’s  function  for  a  source 
of  strength  of  4x 


(7) 


The  gradient  of  the  velocity  potential  function  is  then  forced  to  satisfy  the  boundary  con¬ 
dition  of  (3)  at  a  number  of  control  points  on  the  body.  In  the  current  formulation,  these 
control  points  are  the  panel  centroids,  Figure  1,  and  from  equation  (6),  the  gradient  is  given 
by 


3G(x,xp 

dn 


dT(xj) 


(8) 


The  integrals  in  equations  (6)  and  (8)  can  be  evaluated  exactly,  and  these  analytical  so¬ 
lutions  are  used  when  the  panel- to-control  point  separation  is  relatively  small.  As  the 
separation  increases,  an  expansion  can  be  used  for  the  integral  evaluation  with  little  error. 
The  calculation  of  these  integrals  is  discussed  in  more  detail  in  Appendix  B.  For  a  particular 
set  of  control  points,  the  evaluation  of  the  integrals  in  equations  (6)  and  (8)  leads  to  the 
system 


W  =  [H]W 


i4»  «.•  v, 


{Vn>  =  {I}M  (10) 

with  (/T)  a  coefficient  matrix  relating  source  strength  to  control  point  potential  and  \L) 
a  coefficient  matrix  relating  source  strength  to  control  point  fluid  normal  velocities.  Re¬ 
lations  (9)  and  (10)  define  completely  the  fluid  flow  field  given  the  distribution  of  source 
strengths  a. 

The  fluid  velocities  normal  to  the  panels  at  the  control  points  can  be  related  via  equa¬ 
tion  (3)  to  the  local  nodal  point  velocities  in  the  finite  element  model,  using  the  four  nodes 
which  define  a  quadrilateral  panel,  Figure  1.  For  the  complete  system,  this  becomes 

w.) = (i» 

In  the  current  system,  the  transformation  matrix  [71]  contains  factored  direction  cosines 
between  the  panel  normal  and  the  global  axis  system,  with  each  node  in  the  panel  con¬ 
tributing  1/4  of  the  total.  Equations  (10)  and  (11)  can  be  equated  and  rearranged  to  define 
the  panel  source  strengths  as 

M  =  [L1-'(T]^  (12) 

Substituting  this  expression  into  equation  (9),  we  obtain  an  expression  for  the  velocity 
potential  at  a  control  point  in  terms  of  the  surface  nodal  point  velocities 

{*>  =  (13) 

An  alternative  formulation  for  the  potential  $  can  be  obtained  by  using  a  distribution 
of  dipoles  of  constant  strength  on  a  panel.  The  solution  for  the  dipole  strengths  then  im¬ 
mediately  provides  the  potential  on  the  interface  surface.  In  this  case,  the  integral  equation 
to  be  satisfied  is 

^G(X)X')^1  ^(x')  =  Jf  ,(x>)^  dr(x')  (14) 

in  which  the  derivatives  of  the  velocity  potential  function  are  known.  As  in  the  source 
formulation,  this  can  be  cast  in  discrete  form  as 

(15) 

Solving  for  the  control  point  potential  ind  substituting  for  the  surface  velocities,  we  obtain 

<«  =  m 


Comparing  this  result  with  equation  (13),  it  is  evident  that  the  only  difference  in  the 
two  formulations  is  a  reversal  of  the  order  of  the  coefficient  matrices  in  their  product. 
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Because  both  approaches  represent  fundamentally  exact  solutions,  both  should  converge 
to  the  same  solution  in  the  limit  of  infinite  discretization.  It  is  interesting  to  note  that 
for  certain  structures,  the  added  mass  matrix  obtained  with  a  dipole  formulation  appears 
to  be  significantly  more  accurate  than  that  obtained  from  the  source  formulation.  Similar 
improvements  have  been  discussed  in  the  application  of  the  panel  methods  to  the  study  of 
potential  flow  about  thin  sections.6 

The  pressure  field  in  the  fluid  of  density  p  can  be  determined  from  the  velocity  potential 
via  the  Bernoulli  equation 

p  =  -f>Yt  (1?) 
where  again  the  substantive  derivative  is  used.  For  example,  substitution  of  equation  (13) 
into  (17)  gives  the  panel  control  point  pressures  for  the  source  formulation 


{p>  =  -/>[ff][L]"l[T] 


dt* 


The  pressure  is  integrated  over  the  structure  surface,  in  this  case  by  a  simple  multiplication 
by  the  panel  area  to  give  the  control  point  forces.  These  forces,  which  act  in  the  normal 
direction  to  the  panel,  are  separated  into  global  axis  components  and  distributed  to  the 
appropriate  panel  nodes  via  a  transpose  of  the  transformation  matrix  [T].  The  complete 
expression  for  the  nodal  fluid  force  vector  is  then  given  by 


TO  =  -|riT[x][ff)ii]- 


in  which  [A]  is  a  diagonal  matrix  of  panel  areas.  Neglecting  rigid  body  motion  of  the 
structure  in  the  pressure  and  surface  displacement  derivatives,  and  nonlinear  terms  from 
the  time  derivatives,  (see  Appendix  A)  equation  (19)  reduces  to 

{F}  =  ~[Ma}6  (20) 

with  [Ma]  a  mass  matrix  in  the  form 

[AfA}  =  (r]T[A][^][i]-1(r]  (21) 

for  the  source  formulation,  and 

[Ma)  =  -lT]T[A](L]-1[if][T)  (22) 

for  the  dipole  formulation.  The  terms  in  [Ma]  relate  inertia  forces  at  a  surface  node  in 
a  given  global  direction  to  unit  accelerations  in  a  global  direction  at  a  node.  This  mass 
matrix  can  be  combined  appropriately  with  the  structural  mass  matrix  in  equation  (1) 


If  rigid  body  motion  is  retained  in  the  structural  displacements,  additional  terms  pro¬ 
portional  to  the  surface  velocity  and  displacement  are  present  in  the  nodal  force  vector  (see 
Appendix  A).  The  influence  of  these  terms  on  the  structural  response  does  not  appear  to 
have  been  investigated,  and  is  an  aspect  of  this  approach  worthy  of  further  study.  Their 
influence  is  not  considered  further  in  this  report  since  the  added  damping  terms  cannot  cur¬ 
rently  be  accounted  for  in  the  equilibrium  equation  solution  system,  and  the  added  stiffness 
terms  are  of  higher  order  in  the  displacements  and  are  probably  negligible. 

The  calculation  of  an  added  fluid  mass  matrix  for  a  floating  structure  must  involve  the 
definition  of  the  free  surface.  The  technique  of  imaging  is  used  in  the  surface  panel  method 
to  define  two  forms  of  the  free  surface  boundary  condition  as  well  as  to  define  planes  of 
structural  symmetry.  The  symmetry  option  can  be  used  to  decrease  the  size  of  certain 
problems  where  that  condition  is  applicable. 

For  the  free  surface,  a  pressure  relief  boundary  condition  is  obtained  by  imaging  a  panel 
source  distribution  of  opposite  sign;  a  zero  normal  velocity  surface  condition  is  obtained  by 
imaging  a  distribution  of  the  same  sign.  Equation  (7)  becomes 


for  the  two  free  surface  conditions  available.  In  practice,  the  pressure  relief  condition  pro¬ 
vides  the  correct  approximation  for  most  applications  including  all  those  discussed  in  this 
report.  The  structural  symmetry  condition  is  represented  by  the  form  (25). 

3  Implementation 

The  calculation  of  the  added  fluid  mass  matrix  for  a  finite  element  model  forms  a  small 
part  of  the  complete  system  required  for  dynamic  analyses.  The  organization  of  a  typical 
solution  system  is  shown  in  Figure  2.  Other  than  the  necessary  compatibility  of  surface 
geometry  and  details  of  the  mass  matrix  format,  both  the  finite  element  and  surface  panel 
method  added  mass  systems  are  relatively  independent  of  the  particular  analysis  program. 
The  objective  of  the  current  work  was  to  develop  a  panel  method  as  an  exact  equivalent 
to  the  existing  finite  element  approach,  which  is  discussed  fully  in  Reference  7.  Aside 
from  a  small  increase  in  program  operational  commands  in  the  surface  panel  method,  the 
differences  in  the  systems  are  largely  transparent  to  the  user.  The  differences  in  the  input 
data  requirements  can  be,  of  course,  very  significant. 


3.1  Program  organization 

The  panel  method  calculation  of  the  added  fluid  mass  is  performed  by  the  ADMASS 
program  suite.  The  suite  comprises  four  sections  which  must  be  run  sequentially.  Sec¬ 
tions  can  be  run  independently,  or  combined  as  subroutines  in  a  complete  system.  A  brief 
description  of  each  section  is  given  here. 

Section  1  (ADMS1)  accesses  the  files  containing  structural  geometry  and  finite  element 
fluid  and  fluid/structure  interface  data,  the  format  of  which  is  described  fully  in  Reference 
5.  These  data  are  used  to  define  the  quadrilateral  surface  panels.  The  necessary  panel 
properties  are  calculated  (see  Appendix  B  and  Reference  8),  checks  are  made  on  control 
point  locations,  and  the  properties  stored.  This  section  also  defines  the  [T]  transformation 
matrix  terms,  although  these  are  not  stored  specifically  in  the  [T]  matrix  form  due  to  the 
highly  sparse  nature  of  that  matrix. 

Section  2  (ADMS2)  reads  the  panel  properties  file  and  forms  the  influence  coefficient 
matrices  [ H ]  and  [L\.  The  integral  influence  terms  are  evaluated  either  by  an  exact  integra¬ 
tion,  a  three  term  mulitpole  expansion,  or  a  single  pole  expansion,  depending  on  the  panel 
and  control  point  separation.  This  routine  is  a  modified  version  of  the  influence  coefficient 
routine  currently  used  in  the  program  suite  developed  to  predict  the  potential  flow  about 
ship  hulls.8  Further  details  are  given  in  Appendix  B. 

Section  3  (ADMS3)  computes  the  inverse  of  the  potential  coefficient  matrix  [L] .  This 
system  uses  a  block  diagonal  Gauss-Seidel  iterative  solution  system  which  provides  solution 
vectors  within  a  prescribed  error  tolerance.  This  method  has  been  found  to  be  quite  efficient 
as  an  in-core  solver.  The  block  size  is  currently  a  user-controllable  feature,  since  this 
parameter  can  affect  the  inverse  calculation  efficiency  and,  to  a  lesser  extent,  the  accurac\ 
of  the  inverse  matrix. 

Section  4  (ADMS4)  recovers  the  coefficient  matrices,  the  panel  properties  and  the  trans¬ 
formation  data  and  computes  the  added  mass  matrix  given  in  equation  (21)  or  (22).  This 
matrix  is  stored  as  a  binary  disc  file  in  the  format  required  by  the  VAST  program  for 
addition  to  the  structural  mass  matrix. 

While  the  finite  element  method  provides  a  symmetric  added  fluid  mass  matrix  due 
to  the  variational  basis  of  that  formulation,  the  panel  method  does  not.  Since  the  VAST 
decomposition  algorithm  assumes  symmetry  in  both  stiffness  and  mass  matrices,  a  post¬ 
processing  symmetrizing  option  for  the  ADMASS  system  has  been  implemented  (ADMSP). 
This  system  can  again  be  run  as  a  stand-alone  program,  or  included  in  the  complete  mass 
matrix  generation  process.  The  current  symmetrizing  algorithm  uses  the  simple  formula 

\MA).vm  =  \(\MA)  +  \MA)T)  (26) 

which  is  equivalent  to  a  least  square  minimization  of  the  error  in  the  off-diagonal  terms.  In 
practical  applications,  the  differences  in  results  between  the  symmetric  and  unsymmetric 
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matrices  are  not  always  significant  since  the  mass  matrices  produced  by  the  surface  panel 
method  are  usually  nearly  symmetric. 

All  disc  files  accessed  and  created  by  the  ADMASS  system  are  binary  format,  with  the 
exception  of  the  nodal  coordinate  and  interface  element  data  files  within  the  VAST  system, 
and  a  panel  property  and  diagnostic  data  file  produced  by  the  ADMS1  routine.  These  files 
are  shown  in  the  schematic  of  Figure  3.  Specific  details  on  the  operation  of  the  ADMASS 
system  are  given  in  Appendix  C. 

3.2  Operational  constraints 

Fundamental  differences  between  the  finite  element  and  surface  panel  method  of  calcu¬ 
lating  added  mass  matrices  impose  certain  operational  constraints  in  the  latter’s  application 
to  finite  element-based  models.  Several  key  requirements  have  been  identified  during  the 
implementation  and  testing  of  the  surface  panel  method,  and  these  are  discussed  below. 

3.2.1  Coordinate  systems 

The  imaging  planes  used  for  free  surface  or  athwartship  symmetry  plane  definition  in 
the  ADMASS  system  are  preset  within  the  program;  however,  the  added  nass  matrix  is 
calculated  in  a  coordinate  system  defined  by  the  geometry  of  the  structure.  For  analyses 
in  which  imaging  planes  are  required,  these  internal  and  geometric  coordinate  systems 
must  coincide  in  order  that  the  mass  matrix  produced  can  be  correctly  combined  with  the 
structural  mass  matrix.  In  the  current  ADMASS  system,  the  waterplane  must  be  the  origin 
of  the  global  Z  coordinate  axis  for  free  surface  imaging,  and  the  athwartship  centerplane 
must  be  the  origin  of  the  global  Y  coordinate  axiB  for  symmetry  imaging,  Figure  4.  There  is 
no  provision  for  longitudinal  symmetry.  Since  this  is  not  a  conventional  coordinate  system 
for  structural  definition  (but  is  effective  for  the  imaging  system),  a  geometry  pre-processor 
shown  as  the  first  stage  in  the  flowchart  of  Figure  3  has  been  developed  for  several  of  the 
typical  structural  coordinate  systems.  Note  that  the  entire  VAST  analysis  must  then  be 
done  with  the  geometry  in  this  new  coordinate  system.  For  analyses  in  which  a  free  surface 
or  symmetry  plane  is  not  required,  no  coordinate  transformations  are  necessary. 

3.2.2  Panel  corner  connectivity 

A  second  operational  constraint  derives  from  the  panel  normal  definition  in  the  surface 
panel  method.  These  are  currently  determined  from  vector  cross-products  of  the  panel 
diagonals,  hence  the  panel  nodal  connectivity  is  critical  to  the  correct  orientation  of  the 
normal.  To  obtain  correct  normal  orientation,  a  panel  must  be  defined  by  a  clockwise  corner 
node  order  looking  onto  the  structure  from  the  fluid.  If  required,  an  option  is  available  in 
the  ADMS1  routine  to  reverse  the  complete  set  of  panel  node  ordering  from  that  read  from 


the  element  interface  data  file.  The  VAST  interface  elements  which  the  ADMASS  system 
is  compatible  with  are  shown  in  Figure  5.  The  correct  nodal  connectivity  and  the  method 
of  edge  refinement  for  the  8-noded  element  are  also  shown  in  the  figure. 

Panel  properties  including  areas,  corner  nodes  and  normal  components  can  be  verified 
in  the  ADMS1.DAT  file  produced  by  the  ADMS1  routine,  and  it  is  suggested  that  this 
section  be  run  as  a  stand-alone  to  obtain  this  file  before  the  subsequent  routines  are  run. 
The  ADMS1  section  is  always  the  least  demanding  of  the  four  sections  in  terms  of  CPU 
time,  and  rarely  takes  more  them  a  few  seconds  to  run. 

3.2.3  Non-unique  interface  elements 

A  further  operational  problem  arises  when  an  interface  element  is  connected  to  two 
fluid  finite  elements,  one  on  either  side,  in  the  finite  element  model.  This  is  equivalent 
to  a  zero  thickness  assumption,  and  is  apparently  adequate  for  added  mass  predictions 
by  finite  element  method  for  relatively  thin  structures  such  as  propeller  blades.  For  the 
surface  panel  method,  this  modelling  is  not  allowable,  and  each  surface  must  be  uniquely 
panelled.  To  this  end,  the  ADMS1  routine  offers  several  options  to  create  the  second  surface, 
depending  on  the  type  of  elements  used  in  the  structural  modelling.  In  some  cases,  the  added 
mass  matrix  computed  from  the  two  surface  expansion  must  be  condensed  to  eliminate  the 
second  surface.  This  operation  has  been  included  in  the  same  post-processing  system  as  the 
symmetrizing  option.  For  further  details  of  the  program  options,  see  Appendix  C. 

For  cases  in  which  a  second  surface  cannot  be  defined  with  these  options,  the  ADMASS 
system  will  not  provide  an  added  fluid  mass  matrix  directly.  A  good  approximation  may  be 
obtained  by  simply  multiplying  the  fluid  mass  matrix  determined  with  fluid  on  one  side  by 
a  factor  of  two  in  cases  where  the  surface  normals  on  either  side  are  approximately  collinear. 

4  Application  and  Evaluation 

Subject  to  the  constraints  discussed  in  the  previous  section,  the  surface  panel  method 
provides  an  alternative  capability  for  the  calculation  of  added  fluid  mass  matrices  for  a 
wide  range  of  structural  dynamics  problems.  It  remains  to  evaluate  the  performance  of 
the  system  in  terms  of  the  accuracy  and  efficiency  of  the  method  relative  to  the  existing 
finite  element  system.  There  are  a  number  of  possible  evaluation  procedures;  the  simplest 
is  to  use  finite  element  models  of  structures  for  which  exact  analytical  results  are  available, 
such  as  cylinders  and  spheres,  and  to  compare  the  exact  results  to  the  rigid  body  mass 
predictions  (heave,  sway  and  surge)  from  the  panel  method.  These  rigid  body  terms  are 
obtained  from  a  summation  of  appropriate  terms  in  the  fluid  mass  matrix.  In  this  study,  a 
half-cylinder  has  been  used  for  this  test.  These  results  also  provide  a  good  indication  of  the 
convergence  characteristics  of  the  method  as  a  function  of  the  discretization  of  the  structure 
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surface.  As  a  second  example  of  rigid  body  mass  comparison,  the  heave  and  sway  added 
fluid  mass  predictions  for  two  ship  hulls  are  compared  to  the  predictions  of  strip  theory 
using  the  Lewis  form  method.  For  the  heave  added  mass,  a  comparison  with  a  line  dipole 
distribution  method  is  also  given. 

Since  the  above  comparisons  involve  only  an  integral  of  the  added  fluid  mass,  they  reveal 
nothing  about  the  distribution  of  the  mass  and  its  effect  on  the  dynamic  response  solution. 
To  better  evaluate  the  characteristics  of  the  complete  mass  matrix,  a  number  of  natural 
frequency  analyses  have  been  performed,  utilizing  first  a  fluid  mass  matrix  generated  via  the 
finite  element  formulation,  and  secondly  a  number  of  fluid  mass  matrices  from  the  surface 
panel  method.  In  turn,  these  results  are  compared  to  experimentally  determined  natural 
frequencies.  The  examples  considered  for  these  comparisons  are  a  floating  half-cylinder  and 
two  different  propeller  blades. 

4.1  Half-cylinder  tests 

The  sectional  added  fluid  mass  in  heave  of  a  floating  half-cylinder  of  radius  r  of  infinite 
length  is 

Mk  =  ipxr2  (27) 

and  the  sectional  sway  added  mass  is 

Mlw  =  - pr 2  (28) 

IT 

The  surge  mass  for  the  half-cylinder  considered,  with  a  length/radius  ratio  of  11.4,  can  be 
approximated  adequately  with  the  formula  for  ellipsoids9 

M,e  =  ^kpnr*  (29) 

with  k  =  0.20.  For  r  =4.2  and  L  =96,  the  rigid  body  predictions  from  the  ADMASS  system 
should  converge  to  the  values 

Mh  =  2655  M,w  =  1076  M„  =  61  (30) 

The  added  mass  predictions  from  a  series  of  tests  with  various  circumferential  disretiza- 
tions  of  this  cylinder  are  shown  non-dimensionalized  against  the  above  values  in  Figure  6. 
A  10  panel  longitudinal  discretization  of  the  model  was  used  to  obtain  these  results.  Note 
that  for  these  models  with  fluid  on  one  side  of  the  structure  only,  the  same  results  are 
obtained  from  both  the  source  and  dipole  formulations  of  the  surface  panel  method. 

It  is  evident  that  the  convergence  characteristics  of  the  mass  terms  are  different,  with 
surge  converging  very  quickly,  and  sway  the  most  slowly.  Sway  appears  to  converge  mono- 
tonically  from  above,  while  the  heave  predictions  initially  increase,  then  decrease  through 
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the  theoretical  value.  The  errors  arise  from  two  main  causes:  the  geometric  approximation 
involved  in  a  planar  panel  method  which  results  in  a  poor  representation  of  the  true  surface 
area  of  curved  surfaces,  and  the  fact  that  the  control  points  are  not  on  the  true  surface. 
Similar  convergence  characteristics  were  noted  in  Reference  3  for  a  two-dimensional  source 
distribution  method.  It  is  stated  there  that  the  pressure  terms  are  overestimated  in  the  sur¬ 
face  panel  approach,  causing  the  positive  errors  in  the  estimates  of  sway  and  heave  mass. 
This  may  be  the  case  in  the  ADMASS  system  as  well,  although  the  actual  values  of  those 
terms  have  not  been  investigated.  Based  on  these  results,  the  convergence  of  the  panel 
method  appears  to  be  quite  acceptable.  Tests  with  varying  numbers  of  axial  panels  and  the 
same  circumferential  discretizations  as  Figure  6  gave  very  similar  values.  This  agreement 
correctly  indicates  the  two-dimensional  nature  of  the  rigid  body  mass  distribution  for  a  long 
slender  body. 

The  variety  of  discretization  levels  used  in  these  and  other  tests  have  been  used  to 
construct  the  minimum  CPU  demand  curve  for  the  ADMASS  system,  presented  in  Figure 
7.  These  times  are  for  DREA’s  DEC  2060  machine  and  will  vary  a  small  amount  depending 
on  the  operating  load.  The  actual  run-time  will  also  depend  on  the  convergence  performance 
of  the  matrix  inversion  routine,  which  is  monitored  with  an  iteration  counter  in  the  ADMS3 
routine.  The  curve  in  Figure  7  is  based  on  a  maximum  iteration  count  of  2,  which  is  almost 
always  obtained  for  models  with  fluid  on  one  side  only.  For  thin  submerged  structures,  the 
iteration  count  can  be  higher,  and  the  CPU  demand  will  increase  accordingly. 

The  system  has  also  been  evaluated  on  a  VAX  11-750  machine  where  CPU  demand  will 
average  2.2  -  2.5  times  that  shown  in  Figure  7.  These  CPU  times  do  not  include  any  post¬ 
processing  time,  as  neither  the  symmetric  nor  condensed  form  of  the  added  mass  matrix 
is  required  for  rigid  body  estimates.  The  current  system  utilizes  several  in-core  solution 
algorithms  for  efficiency,  thus  imposing  a  limit  on  the  matrix  size  which  can  be  processed. 
For  the  DEC  2060,  this  limit  is  approximately  650  panels. 

The  rigid  body  estimates  for  any  structure  are  obtained  from  the  addition  of  all  mass 
components  in  a  particular  degree  of  freedom  in  the  added  mass  matrix.  For  the  cylinder 
model  for  example,  heave  mass  is  obtained  from  the  addition  of  every  third  element  in  every 
third  row  of  the  matrix.  Sway  and  surge  masses  can  be  calculated  similarly. 

4.2  Ship  hull  added  mass  comparison 

The  ADMASS  system  has  been  used  to  compute  rigid  body  added  mass  matrices  for 
two  frigate  hulls  for  which  comparisons  can  be  made  with  the  Lewis  form  method10  and  a 
line  dipole  distribution  method.11  The  models  used  in  these  analyses  were  generated  with 
a  modified  version  of  the  hull  surface  generating  routines  developed  for  the  prediction  of 
potential  flow  around  ship  hulls.12,13  This  system  produces  half-hull  models,  which  can  be 
used  with  the  centerplane  symmetry  option  in  the  ADMASS  system.  To  test  this  option, 
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Table  1:  Added  fluid  mass  comparison  for  the  DDH265  and  CPF  hulls. 


the  half-hull  models  were  mirrored  to  obtain  a  full  hull  model,  and  results  compared  for  the 
two  approaches.  The  inverted  full  hull  surfaces  panelled  with  200  elements  are  shown  in 
Figures  8  and  9.  The  ADMASS  results  for  both  the  DDH265  and  Canadian  Patrol  Frigate 
are  presented  in  Table  1,  as  are  the  estimates  of  added  mass  from  the  other  methods.  A 
consistent  draft  has  been  maintained  in  the  various  added  mass  evaluations.  The  Lewis 
form  results  were  generated  using  the  sectional  area  and  half-beam  and  draft  values  for  21 
stations  on  each  of  the  hulls.  Those  predictions  do  not  include  a  modal  J  reduction  factor 
to  account  for  three-dimensional  effects. 

These  results  indicate  that  the  symmetry  option  correctly  produces  one-half  of  the  total 
mass  for  the  heave  and  surge  directions.  The  sway  results  are  incorrect  for  the  half- models; 
however,  the  full  model  results  are  considered  accurate.  The  cause  of  the  discrepancy  in 
the  half-model  results  is  under  investigation.  On  the  basis  of  this  comparison,  the  surface 
panel  method  appears  to  overpredict  the  rigid  body  heave  and  sway  mass  with  convergence 
from  above,  consistent  with  the  results  of  the  floating  cylinder  tests.  The  agreement  in  the 
results  is  quite  encouraging,  since  these  models  represent  relatively  coarse  grids.  In  the  case 
of  sway,  the  Lewis  form  methods  are  not  considered  to  be  particularly  accurate;  the  panel 
method  results  for  the  full  model  are  probably  closer  to  the  actual  values. 

4.3  Natural  modes  of  a  floating  cylinder 

To  provide  a  more  comprehensive  evaluation  of  the  characteristics  of  the  added  mass 
matrix  produced  by  the  panel  method,  the  fluid  mass  matrices  from  the  ADMASS  and 
finite  element  formulation  have  been  used  in  the  natural  frequency  analysis  of  a  floating 
half-cylinder.  This  finite  element  model  also  represents  a  physical  model  for  which  natural 
frequencies  have  been  determined  experimentally14  for  various  drafts.  The  comparisons  of 
the  frequencies  of  this  model  for  one  draft  are  presented  in  Table  2.  The  finite  element 


Vibration  Mode 

2  Layer  F.E.M. 

Panel  Method 

(N  =  node) 

Hz 

Hz 

2N  Vertical 

7.8 

7.20 

7.8 

3N  Vertical 

20.8 

19.8 

21.3 

4N  Vertical 

40.0 

38.4 

40.8 

5N  Vertical 

64.1 

62.6 

65.8 

2N  Horizontal 

7.8 

7.20 

7.8 

3N  Horizontal 

30.8 

13.5 

16.0 

4N  Horizontal 

36.6 

25.4 

27.8 

5N  Horizontal 

52.7 

41.1 

43.8 

IN  Torsional 

1.6 

5.5 

2N  Torsional 

36.6 

41.5 

3N  Torsional 

30.8 

16.0 

4N  Torsional 

52.7 

27.8 

5N  Torsional 

- 

41.1 

43.8 

Surge 

0.145 

0.166 

Sway 

4.37 

5.55 

Heave 

26.8 

33.7 

Table  2:  Comparison  of  natural  mode  predictions  for  a  floating  cylinder. 

fluid  model  used  two  layers  of  20-node  fluid  elements,  and  approximately  5.5  hours  of  CPU 
time  were  required  to  generate  the  added  mass  matrix.  A  three  layer  fluid  model  was  also 
tested,  which  gave  results  very  similar  to  those  of  the  two  layer  model.  This  is  consistent 
with  finite  element  fluid  modelling,  in  which  two  layers  of  quadratic  elements  will  often 
adequately  represent  the  infinite  fluid  domain.  A  three  layer  fluid  model  for  this  structure 
does  not  represent  a  particularly  tractable  problem  due  to  the  large  CPU  time  required  for 
the  generation  of  the  added  mass  matrix. 

The  rigid  body  added  mass  predictions  for  the  two  methods  are  also  shown  at  the 
bottom  of  Table  2.  The  ADMASS  system  gives  significantly  higher  rigid  body  added  mass 
estimates  than  the  finite  element  method;  this  is  shown  to  be  a  consistent  result  in  other 
comparisons  of  this  nature.  The  ADMASS  CPU  time  for  the  fluid  mass  matrix  generation 
for  this  168  panel  model  is  about  4.5  minutes.  The  vertical  natural  frequency  predictions 
utilizing  either  the  finite  element  or  surface  panel  method  are  in  good  agreement  with 
the  experimental  results.  The  larger  mass  from  the  panel  method  indicated  in  the  rigid 
body  comparisons  appears  to  result  in  a  slightly  low  prediction  of  these  vertical  modes. 
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In  the  horizontal  and  torsional  modes,  the  panel  method  mass  results  appear  to  provide 
better  overall  accuracy  than  those  obtained  with  the  finite  element  mass  matrix.  The  exact 
cause  of  the  poor  performance  of  the  finite  element  model  in  the  prediction  of  these  modes 
can  probably  be  attributed  to  the  complex  coupling  which  occurs  between  the  horizontal 
and  torsional  displacements,  combined  with  the  inability  of  the  shell  element  used  in  the 
structural  modelling  to  properly  represent  the  shear  stiffness  distribution  in  this  thin-walled 
structure.  The  differences  between  the  results  for  t>*e  panel  method  mass  matrix  and  the 
finite  element  mass  matrix  is  a  result  of  the  sensitivity  of  this  structure  to  the  fluid  mass 
distribution,  which  is  a  dominant  component  in  the  overall  mass  matrix. 

It  should  be  noted  that,  for  this  case  and  any  other  that  utilizes  the  8-noded  interface 
element,  the  mass  matrix  produced  by  the  ADMASS  system  is  significantly  smaller  than 
that  produced  by  the  finite  element  method.  This  disparity  results  from  the  manner  in 
which  panels  can  be  defined  from  elements  which  include  mid-side  nodes,  Figure  5.  The 
default  option  takes  only  the  four  corner  nodes  of  this  element,  thus  producing  a  mass 
component  for  only  half  of  the  interface  node  .  An  option  to  refine  either  edge  of  the  basic 
grid  using  two  of  the  mid-side  nodes  is  available  within  the  ADMS1  routine.  This  option 
then  provides  a  mass  component  for  6  of  the  8  nodes  in  each  interface  element.  Currently 
there  is  no  provision  for  a  panelling  which  provides  a  mass  component  to  all  8  nodes  of 
the  interface  element,  since  this  would  require  the  creation  of  a  new  node  within  the  panel. 
The  contribution  of  that  node  would  then  have  to  be  condensed  out  of  the  final  added  mass 
matrix  before  it  could  be  combined  with  the  structural  mass  matrix.  Mass  components  for 
all  interface  nodes  are  provided  when  the  4-node  interface  element  is  used.  For  typical  finite 
element  grids,  it  is  unlikely  that  a  refinement  further  than  that  available  will  be  required 
for  accurate  assessment  of  the  added  fluid  mass.  In  the  above  example,  the  default  option 
was  used. 

4.4  Natural  modes  of  propellers 

In  a  similar  manner  to  the  above  example,  the  added  mass  matrices  produced  by  the 
finite  element  and  surface  panel  methods  have  been  used  in  natural  frequency  analyses  of 
two  propeller  blades.  These  models  again  provide  a  good  comparison  basis  because  of  the 
large  percentage  of  the  total  mass  which  stems  from  the  fluid  effects.  Changes  in  this  mass 
distribution  thus  have  fairly  evident  effects  on  the  natural  frequency  predictions.  The  first 
analysis  is  of  a  destroyer  blade  which  has  little  skew  and  presents  no  particular  modelling 
difficulties.  The  second  is  of  an  experimental,  model-sized  blade  having  very  high  skew. 
The  finite  element  models  of  these  blades  are  shown  in  Figures  10  and  11.  A  mid-surface 
view  of  the  3-dimensional  fluid  model  is  also  shown  surrounding  the  blades.  The  structural 
and  fluid  models  were  created  with  a  propeller  modelling  program.15 


4.4.1  Destroyer  Blade 

The  destroyer  blade  served  as  a  test  case  for  the  comparison  of  several  different  modelling 
options  and  discretization  levels  and  for  the  evaluation  of  post-processing  features  within 
the  ADMASS  system.  Two  structural  modelling  options  can  be  used  in  the  current  propeller 
modelling  system:  a  shell  element  defined  by  8  displacement  nodes  on  one  surface  and  8 
geometric  nodes  on  the  opposite  surface,  or  a  20-node  solid  element.8  Both  of  these  models 
utilize  the  8-node  fluid/structure  interface  element,  but  only  one  surface  of  the  blade  is 
used  as  the  interface,  with  fluid  elements  attached  to  either  side  of  each  interface  element. 
This  thin  body  assumption  is  not  compatible  with  the  surface  panel  method,  but  for  either 
modelling  option,  it  is  possible  to  define  the  required  second  surface  based  on  the  geometry 
of  the  interface  elements  and  a  knowledge  of  the  node  numbering  characteristics  of  the 
model  generation  routine. 

In  the  case  of  the  20-node  solid  model,  the  resulting  fluid  mass  matrix  from  the  panel 
method  is  larger  than  that  produced  by  the  finite  element  method  with  the  thin  body 
assumption.  That  matrix  can  be  combined  directly  with  the  structural  mass  matrix,  since 
all  nodes  involved  are  displacement  nodes.  For  the  16-node  shell  modelling  option,  mass 
components  on  the  second  surface  are  associated  with  geometric  nodes,  which  are  not 
included  in  the  overall  analysis.  These  mass  components  must  therefore  be  condensed 
out  of  the  total  added  fluid  mass  matrix  by  combination  with  the  associated  displacement 
node  terms.  A  post-processing  system  to  automatically  perform  this  reduction  has  been 
developed  and  must  be  used  before  the  added  fluid  mass  matrix  can  be  combined  with 
the  structural  mass  matrix.  The  reduction  and  symmetrizing  options  are  available  in  the 
routine  ADMSP  described  further  in  Appendix  C. 

The  numerous  comparisons  for  this  propeller  blade  are  summarized  in  Tables  3,  4,  and 
5.  The  results  for  the  solid  models  are  presented  in  Table  3,  and  those  for  the  shell  elements 
in  Table  4.  Results  for  the  source  and  dipole  formulations  are  also  presented  in  those  tables. 
The  fluid  model  for  the  finite  element  added  mass  matrix  for  both  models  consisted  of  two 
layers  of  20-node  fluid  elements,  as  indicated  in  Figure  10.  A  three  layer  model  was  also 
tested  with  one  analysis  of  the  shell  model;  those  results  were  virtually  identical  to  the 
predictions  using  the  two-layer  model,  hence  the  latter  is  considered  an  adequate  basis  for 
comparison  of  accuracy  and  CPU  time.  Experimentally  determined  frequencies16  are  also 
shown  in  Tables  3  and  4. 

It  is  evident  that  several  of  the  measured  frequency  values  are  not  predicted  with  the 
numerical  models.  The  fourth  and  eighth  modes  are  missed  in  all  of  the  solid  models, 
and  the  shell  models  miss  the  fourth  mode  in  several  cases.  These  defficiencies  point  out 
the  sensitivity  of  the  eigenvalue  analysis  to  both  the  stiffness  representation  and  the  mass 
distribution  for  this  structure.  As  can  be  seen  in  Table  4,  the  analysis  of  this  blade  using 
a  more  detailed  grid  of  80  shell  elements  (160  panels)  gave  a  prediction  for  all  of  the 


Mode 


E: 


F.E.  Mass 

Source  panel 

Dipole  panel  | 

50 

50 

100  rad1 

100  rad 

100  chd 

50 

100  rad 

33 

30 

30 

30 

31 

35 

35 

76 

76 

74 

75 

80 

81 

77 

101 

111 

93 

97 

101 

102 

99 

156 

_ 

151 

158 

172 

163 

151 

172 

185 

171 

175 

182 

183 

171 

233 

223 

213 

226 

245 

227 

217 

289 

293 

282 

296 

309 

289 

271 

1160 

23 

100 

100 

100 

23 

100 

1  33  30  3C 

2  76  76  74 

3  101  111  93 

4  - 

5  156  -  15 

6  172  185  17 

7  233  223  21 

8  - 

9  289  293  28! 

CPU  sec*  1160  I  23  I  10< 

[symmetric  form 
lss  matrix  generation  time  only 


:  Comparison  of  natural  mode  predictions  (in  Hz)  for  the  solid  finite  element  model 
estroyer  propeller  blade;  'rad’  and  ‘chd’  refer  to  radial  and  chordwise  respectively. 
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F.E.  Mass 
50 

Source  pam 
50  100  rad 

81 

160 

50 

Dipole  panel 

100  rad  |  100  chd 

Exp1 

32 

29 

30 

31 

34 

34 

35 

36 

73 

77 

73 

74 

80 

76 

84 

79 

90 

95 

92 

91 

97 

94 

98 

96 

128 

- 

133 

127 

- 

136 

- 

131 

152 

149 

- 

160 

151 

158 

153 

150 

178 

185 

167 

198 

172 

179 

174 

178 

224 

218 

206 

209 

184 

- 

186 

206 

- 

255 

228 

251 

253 

228 

251 

240 

300 

327 

290 

293 

311 

286 

310 

280 

1160 

30 

130 

380 

30 

130 

130 

1  32  29  30 

2  73  77  73 

3  90  95  92 

4  128  -  133 

5  152  149 

6  178  185  167 

7  224  218  206 

8  -  255  228 

9  300  327  29C 

CPU  sec1  1  1160  30  13C 

Mass  matrix  generation  time  only 


:  Comparison  of  natural  mode  predictions  (in  Hz)  for  the  shell  finite  element  model 
estroyer  propeller  blade. 
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experimentally  determined  modes.  In  general,  the  solution  utilizing  the  finite  element- 
based  fluid  mass  matrix  gave  very  accurate  frequency  estimates  for  these  models.  The  best 
predictions  were  obtained  when  shell  elements  were  used  to  model  the  propeller  blade. 

The  results  of  six  separate  analyses  to  test  the  ADMASS  fluid  mass  matrix  with  the 
solid  model  are  presented  in  Table  3.  First,  several  analyses  using  the  source  formulation 
were  performed,  starting  with  the  unrefined  finite  element  interface  grid,  which  produces  50 
panels,  25  on  each  face  of  the  propeller.  This  mass  matrix  was  used  both  in  the  unsymmet- 
ric  and  symmetric  form  for  the  source  formulation.  The  results,  which  show  inconsistent 
accuracy  and  missed  modes,  were  the  same  for  either  form  of  the  matrix.  Two  panel  refine¬ 
ment  options,  a  radial  and  chordwise  refinement  of  order  two,  were  also  tested  with  the  solid 
element  model.  It  was  anticipated  that  a  radial  refinement  would  be  the  more  appropriate 
for  a  typical  propeller  blade  because  of  the  radially  biased  aspect  ratio  of  the  structure. 
The  radially  refined  results  do  indicate  slightly  better  agreement  than  those  obtained  with 
the  chordwise  refinement.  The  radially  refined  model  was  also  tested  in  the  unsymmetric 
and  symmetric  form.  Again,  no  significant  change  in  overall  accuracy  was  obtained  with 
the  symmetric  form,  indicating  that  this  is  not  a  major  disadvantage  of  the  surface  panel 
method.  For  consistency,  all  further  runs  utilized  a  symmetric  form  of  the  fluid  mass  matrix. 

The  frequency  predictions  for  the  fundamental  mode  using  the  source  formulation  form 
of  the  added  mass  matrix  are  somewhat  low,  indicating  an  overprediction  of  the  added  mass 
effect.  This  overestimate  is  evident  in  the  rigid  body  added  mass  estimates  produced  by  the 
source  formulation,  which  are  presented  in  Table  5.  The  mass  estimates  tend  to  decrease 
with  increasing  discretization,  consistent  with  the  convergence  characteristics  found  in  the 
half-cylinder  tests,  but  are  in  general  too  high.  The  added  mass  estimates  from  the  dipole 
formulation  are  in  much  better  agreement  with  the  finite  element  fluid  mass  estimates.  Since 
the  finite  element  method  itself  will  overpredict  the  mass  effect  for  a  submerged  body,  these 
results  indicate  that  the  dipole  formulation  is  the  more  accurate  of  the  two  panel  method 
approaches  for  this  structure.  Two  of  the  natural  frequency  analyses  were  repeated  using 
the  dipole  form  of  the  added  mass  matrix:  the  original  grid  and  the  radially  refined  grid. 
Those  results,  given  in  Table  3,  provide  improved  accuracy  for  the  lower  vibration  modes 
of  this  blade. 

As  one  test  to  ensure  that  the  refinement  algorithm  was  functioning  correctly,  the  rigid 
body  fluid  mass  estimates  from  the  radial  refinement  of  the  50  panels  to  100  panels  was 
compared  to  a  model  which  was  generated  with  50  elements  (100  panels).  These  results 
are  shown  in  Table  5,  and  indicate  that  similar,  but  not  exact  agreement  is  obtained.  The 
small  disparity  is  due  to  the  slightly  different  panel  grids  which  are  produced  by  the  two 
methods. 

Table  4  presents  the  results  of  a  similar  series  of  analyses  performed  with  the  shell 
element  models.  The  frequency  predictions  of  the  lower  modes  are  again  low  due  to  the 
overprediction  of  the  mass  effect  for  this  blade  by  the  source  formulation.  The  finite  element- 


Dipole  form  mass  (Kg) 

Surge 

Sway 

Heave 

1020 

337 

20.7 

1044 

391 

23.0 

952 

360 

21.8 

1019 

374 

23.9 

994 

362 

23.4 

890 

327 

22.1 

Model 


F.E.  Mass 
50  panels 
100  chord 
100  radial 
100  panels 
160  panels 


Table  5:  Comparison  of  rigid  body  mass  predictions  from  the  finite  element  and  panel 
method  mass  matrices  for  the  destroyer  propeller  blade. 

based  fluid  mass  matrix  and  that  from  the  dipole  formulation  of  the  panel  method  provide 
about  the  same  overall  accuracy,  although  modes  are  missed  sporadically  in  both  cases. 
Prediction  of  all  experimentally  measured  modes  was  obtained  with  the  more  detailed  shell 
model,  at  the  expense  of  increased  CPU  time.  These  results  have  all  utilized  a  symmetric 
form  of  the  fluid  mass  matrix.  The  first  four  modal  displacement  contours  for  this  propeller 
are  shown  in  Figures  12  -  15  for  the  two  fluid  mass  matrix  options.  The  comparisons  are 
based  on  the  results  from  the  shell  finite  element  model  for  the  cftie  of  the  two-layer  finite 
element  fluid  model  and  the  100  panel  radially  refined  model  using  the  dipole  formulation 
(columns  2  and  7  in  Table  4).  The  agreement  between  the  results  for  the  two  mass  matrices 
is  very  good  for  the  lower  modes.  The  frequencies  of  the  higher  modes  are  also  predicted 
well  with  both  forms  of  the  fluid  mass  matrix,  but  the  modal  displacement  contours  begin  to 
show  differences.  No  experimentally  determined  modal  displacement  contours  are  available 
for  comparison. 

As  discussed,  the  use  of  the  symmetrized  mass  matrix  has  little  effect  on  the  modal 
predictions,  and  it  is  debatable  whether  this  option  is  required,  although  it  is  not  a  de¬ 
manding  one  in  terms  of  CPU  time.  As  the  discretization  level  increases,  the  effect  of  the 
nonsymmetric  form  should  decrease,  although  this  has  not  been  investigated  in  detail. 

The  CPU  time  requirements  presented  at  the  bottom  of  Tables  3  and  4  are  indicative 
of  the  increased  speed  with  which  the  surface  panel  method  can  produce  the  added  fluid 
mass  matrix,  even  when  post- processing  options  are  required.  This  post-processing  can  be 
a  significant  portion  (~  20  percent)  of  the  total  fluid  mass  matrix  generation  time.  For 
the  propeller  analyses,  the  fluid  mass  matrix  generation  by  finite  element  method  takes 
approximately  50  percent  of  the  total  analysis  CPU  time.  The  reduction  by  a  factor  of 
10  of  that  segment  thus  represents  a  large  saving  in  overall  time  and  cost.  An  even  larger 
reduction  ratio  occurs  for  analyses  on  the  VAX  11-750  machine,  where  processing  times  are 
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typically  2.5  times  longer  for  any  computationally  intensive  problems.  These  time  estimates 
do  not  reflect  any  savings  in  model  generation  time,  which  is  essentially  negligible  for  the 
surface  panel  method,  since  it  could  work  directly  from  the  structural  element  definition  for 
a  submerged  body  such  as  the  propeller.  As  well,  the  accuracy  of  the  added  mass  matrix 
does  not  depend  on  the  fluid  modelling,  which  can  present  difficulties  in  automation  for 
more  complex  geometries.  The  CPU  time  required  for  the  overall  analysis  of  the  detailed 
model  including  the  panel  method  added  mass  matrix  was  still  significantly  lower  that  that 
needed  for  analysis  of  the  smaller  models  with  the  finite  element  fluid  mass  matrix. 

4.4.2  DTNSRDC  P4388  Blade 

As  a  final  comparison,  a  highly  skewed,  model-size  propeller  blade  has  been  analyzed  in 
a  similar  manner  as  described  for  the  destroyer  blade.  This  model  represents  a  somewhat 
more  complex  geometry,  and  a  grid  of  45  shell  elements  (16  node  option)  has  been  used  in 
the  finite  element  model,  Figure  11,  in  place  of  the  25  elements  used  for  the  destroyer  blade. 
The  results  of  these  analyses  are  presented  in  Table  6.  The  unrefined  and  both  chordwise 
and  radial  refinement  results  for  the  source  formulation  are  compared  to  those  utilizing  a 
2-layer  fluid  model  for  added  fluid  mass  matrix  creation,  Figure  11,  as  well  as  to  natural 
frequencies  determined  experimentally  by  holographic  methods.17,18  The  results  of  analyses 
of  the  unrefined  and  radially  refined  models  with  the  fluid  mass  matrix  derived  with  the 
dipole  formulation,  and  an  analysis  using  80  shell  elements  and  a  grid  of  320  panels  are 
also  given  in  the  table.  In  this  case,  the  finite  element  method  predicts  a  frequency  which 
appears  to  have  been  missed  in  the  experimental  measurements. 

The  results  for  the  lower  modes  of  this  blade  are  again  quite  encouraging  for  both 
the  finite  element-based  and  source  formulation  panel  method  fluid  mass  matrices.  The 
predictions  using  the  mass  matrix  based  on  the  dipole  formulation  are  poor,  even  for  the 
lower  modes.  The  higher  mode  predictions  of  all  the  analyses  deviate  from  those  determined 
experimentally;  however,  the  agreement  is  within  10  percent  for  most  modes,  which  is 
reasonable.  The  particularity  poor  performance  of  the  dipole  formulation  here,  in  contrast 
to  the  better  accuracy  it  provided  in  the  previous  example,  is  probably  a  result  of  the  very 
small  thickness  of  this  propeller  blade.  Although  this  appears  to  contradict  the  suggestion 
that  the  dipole  formulation  offers  increased  accuracy  for  thin  sections,  it  is  probable  that 
this  model  is  approaching  the  limits  for  which  any  panel  method  could  be  expected  to 
perform  well.  In  these  cases,  the  finite  element-based  added  mass  matrix  may  provide  more 
consistent  accuracy. 

The  rigid  body  added  mass  estimates  for  this  blade  for  the  various  models  are  presented 
in  Table  7.  As  could  be  anticipated  from  consideration  of  the  frequency  predictions,  the 
dipole  formulation  significantly  underestimates  the  dominant  surge  added  mass  component. 
The  convergence  characteristics  for  this  blade  are  also  inconsistent  with  previous  results, 
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Mode  F.E.  Mass  St 

90  90  180 

1  80  83  8' 

2  237  277  23 

3  498  554  5C 

4  590  652  68 

5  897  1026  91 

6  995  1121  10! 

7  1485  1703  14! 

8  1611  1851  171 

CPU  sec1  5900  115  53 

1  Mass  matrix  generation  time  only 
All  results  based  on  a  shell  model 


Table  6:  Comparison  of  natural  mode  predictions  (in  Hz)  for  the  DTNSRDC  P4388  propeller 
blade. 


F.E.  Mass 

Source  panel  | 

Dipole  panel 

Exp1* 

90 

90 

180  rad 

160 

320  rad 

90 

180  rad 

80 

83 

84 

80 

80 

92 

90 

82 

237 

277 

236 

218 

225 

296 

266 

232 

498 

554 

506 

477 

472 

594 

543 

- 

590 

652 

688 

690 

600 

675 

619 

576 

897 

1026 

919 

857 

849 

1055 

941 

738 

995 

1121 

1090 

1079 

956 

1139 

1017 

930 

1485 

1703 

1497 

1360 

1345 

1765 

1506 

1245 

1611 

1851 

1761 

1633 

1630 

1843 

1630 

1600 

5900 

115 

535 

380 

2621 

115 

540 

Model 

Source  form  mass  (Kg) 
Surge  Sway  Heave 

Jipole  form  mass  (Kg) 
Surge  Sway  Heave 

F.E.  Mass 

0.602 

0.127 

0.093 

0.602 

0.127 

0.093 

90  panels 

0.569 

0.133 

0.117 

0.465 

0.114 

0.097 

160  panels 

0.701 

0.178 

0.117 

0.481 

0.121 

0.095 

180  radial 

0.667 

0.158 

0.118 

0.502 

0.175 

0.126 

180  chord 

0.632 

0.150 

0.114 

0.473 

0.118 

0.094 

320  radial 

0.666 

0.155 

0.123 

0.472 

0.120 

0.090 

Table  7:  Comparison  of  rigid  body  mass  predictions  from  the  finite  element  and  panel 
method  mass  matrices  for  the  P4388  propeller  blade. 
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again  suggesting  that  the  panel  method  is  having  difficulty  with  the  geometry  of  this  model. 
As  well,  the  inaccuracy  of  the  higher  mode  predictions  in  all  the  analyses  indicates  the 
sensitivity  of  the  overall  analysis  to  the  numerical  stiffness  and  mass  representations.  The 
results  for  the  highly  discretized  model  show  improved  accuracy,  but  the  CPU  time  required 
for  that  analysis  is  unacceptably  high. 

Four  of  the  modal  displacement  contours  for  this  blade  are  presented  in  Figures  16  -  19. 
These  comparisons  are  based  on  the  finite  element  fluid  model  and  the  radially  refined  model 
(180  panels)  with  the  source  formulation  fluid  mass  matrix.  There  is  excellent  agreement 
in  the  lowest  two  modes,  but  slight  differences  become  evident  in  the  higher  modes.  This 
again  highlights  the  sensitivity  of  the  analyses  to  the  fluid  mass  distribution. 

The  CPU  times  for  generation  of  the  various  fluid  mass  matrices  used  in  the  comparisons 
for  this  blade  are  shown  at  the  bottom  of  Table  6.  The  times  given  for  the  surface  panel 
method  include  the  condensation  and  symmetrizing  time.  As  was  noted  in  the  previous 
example,  a  saving  of  about  an  order  of  magnitude  in  CPU  time  can  be  obtained  through 
the  use  of  the  surface  panel  source  method  in  place  of  the  two  layer  fluid  models. 

5  Concluding  Remarks 

The  example  applications  have  shown  that  the  surface  panel  method  is  an  effective 
alternative  to  the  current  method  of  finite  element  fluid  modelling  for  the  generation  of  the 
added  fluid  mass  matrix.  In  all  cases,  the  surface  panel  method  provided  a  fluid  mass  matrix 
which  gave  similar  accuracy  for  the  lower  vibration  modes  in  significantly  less  CPU  time 
than  that  required  for  matrix  generation  via  the  finite  element  method.  For  moderately 
thin  submerged  structures  such  as  full  scale  propeller  blades,  the  dipole  formulation  of  the 
surface  panel  method  provides  better  overall  accuracy  than  the  source  formulation.  For 
very  thin  structures,  the  reverse  appears  to  be  true,  and  the  finite  element  method,  in 
conjunction  with  a  thin  body  assumption,  will  probably  provide  more  consistent  accuracy 
than  either  form  of  the  surface  panel  method. 

The  rigid  body  added  mass  estimates  from  the  surface  panel  method  for  the  structures 
with  fluid  on  only  one  side  indicate  that  the  method  will  overpredict  the  added  mass  effect, 
and  converge  with  increasing  discretization  from  above  the  exact  value.  This  results  in  a 
convergence  from  below  for  the  lowest  vibration  modes.  The  overprediction  has  little  effect 
on  the  higher  vibration  modes.  For  moderately  thin,  submerged  structures,  the  convergence 
characteristics  follow  similar  trends,  but  tend  to  be  more  sensitive  to  the  panel  geometry, 
and  may  show  inconsistencies.  It  is  suggested  that  this  dependence  on  interface  geometry 
is  no  more  severe  than  that  for  the  finite  element  method  for  such  structures,  although  that 
sensitivity  has  not  been  investigated  in  detail  in  the  current  work.  For  very  thin  structures, 
these  inconsistencies  can  be  acute,  particularity  for  the  dipole  formulation.  The  use  of  finite 


element-based  added  mass  matrices  will  provide  greater  accuracy  in  these  cases,  since  that 
method,  in  conjunction  with  the  thin  body  assumption,  will  not  be  affected  by  changes  in 
the  structural  thickness. 

A  finite  element  grid  adequate  for  vibration  analyses  of  the  structure  in  air  generally 
provides  a  definition  of  hydrodynamic  panelling  adequate  for  the  prediction  of  the  lower 
vibration  modes  of  the  structure  in  water.  For  accuracy  in  the  higher  mode  predictions, 
refinement  of  the  basic  grid  will  be  required.  This  refinement  may  have  an  adverse  effect 
on  the  convergence  performance  of  the  coefficient  matrix  inversion  algorithm  for  structures 
with  reasonably  complex  geometry.  The  convergence  performance  can  usually  be  improved 
by  increasing  the  user-controlled  block  size  in  that  routine. 

The  lack  of  symmetry  in  the  fluid  added  mass  matrices  generated  with  the  surface  panel 
method  appears  to  have  little  effect  on  the  accuracy  of  the  vibration  mode  predictions  for 
the  structures  analysed.  A  simple  and  non-CPU-intensive  symmetrizing  algorithm  can  be 
used  which  provides  a  slight  increase  in  overall  accuracy. 

For  the  models  considered,  the  CPU  time  ratios  for  the  finite  element  method  versus  the 
panel  method  varied  from  a  minimum  of  approximately  10  in  the  propeller  analyses,  to  a 
maximum  of  about  70  for  the  floating  cylinder.  This  ratio  is  necessarily  very  dependent  on 
the  fluid  model  used  for  the  comparison,  but  in  general,  a  CPU  time  reduction  of  at  least 
an  order  of  magnitude  can  be  expected  using  the  surface  panel  method.  For  large  and/or 
complex  structures  for  which  the  necessity  of  fluid  modelling  has  previously  hampered  the 
analyses,  the  surface  panel  method  may  provide  an  effective  fluid  mass  modelling  capability. 

Several  aspects  of  the  surface  panel  method  for  added  fluid  mass  prediction  could  be 
examined  in  more  detail  or  developed  to  provide  increased  accuracy.  Typically,  higher 
accuracy  will  come  at  the  expense  of  increased  computation  time.  Two  general  areas  are 
suggested  for  further  investigation  and  development: 

1.  the  development  of  more  flexible  methods  of  defining  panels  from  finite  element  grid 
data,  and; 

2.  the  incorporation  of  a  more  general  Green’s  function  in  the  evaluation  of  the  influence 
coefficients. 

Addressing  the  first  point,  a  large  part  of  the  development  of  a  practical  panel  method 
system  has  involved  the  conversion  of  finite  element  grid  specifications  to  effective  hydro¬ 
dynamic  panelling.  Further  effort  could  be  directed  at  providing  more  flexibility  in  the 
panel  definition  in  terms  of  both  refinement  and  enlargement  options,  since  there  is  no 
requirement  that  the  elements  and  panels  be  related  on  a  one-to-one  basis.  Refinement 
beyond  the  level  currently  available  in  the  system  will  create  additional  nodes,  for  which 
the  mass  contribution  will  have  to  be  condensed  out  of  the  final  fluid  mass  matrix.  This 
option  would  allow  the  generation  of  a  fluid  mass  matrix  with  mass  contributions  for  all 
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nodes  of  the  8-node  interface  element.  An  enlargement  option  might  be  convenient  in  cases 
where  a  detailed  model  created  for  stress  analysis  is  used  in  a  vibration  analysis,  where  such 
refinement  may  not  be  required. 

The  second  important  area  of  further  research  involves  the  replacement  of  the  simple  free 
space  Green’s  function  with  a  frequency  dependent  function  which  accounts  for  the  presence 
of  waves  on  the  free  surface.19'30  The  implementation  of  that  option  from  a  programming 
viewpoint  would  be  quite  straightforward,  although  the  more  complex  function  will  require 
a  numerical  integration  algorithm  in  addition  to  the  exact  and  multipole  methods  exploited 
in  the  current  system.  It  is  anticipated  that  this  approach  would  have  adverse  effects  on 
the  computational  efficiency  of  the  panel  method  system,  but  could  increase  the  accuracy 
in  cases  where  free  surface  effects  are  important. 

In  conclusion,  the  inherent  advantages  of  the  surface  panel  method  for  added  mass 
matrix  calculation  can  be  realized  within  a  finite  element-based  analysis  system.  With  the 
exception  of  very  thin  structures,  the  panel  method  should  provide  an  effective  alternative 
to  finite  element  modelling  of  infinite  fluid  domains. 
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Figure  1:  Surface  panel  nomenclature  and  representation. 


Figure  2:  Dynamic  response  analysis  flow  chart  for  structures  requiring  added  fluid  mass 
evaluation. 


Figure  3:  ADMASS  system  flowchart. 
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Figure  4:  Internal  coordinate  system  for  imaging  in  the  ADMASS  system. 
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Figure  5:  Interface  elements,  nodal  connectivity  and  refinement  options  for  the  ADMASS 
system. 
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Figure  6:  Rigid  body  mass  convergence  for  &  half-cylinder  as  a  function  of  circumferential 
discretization. 
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Figure  7:  Minimum  CPU  time  demand  for  mass  matrix  generation  as  a  function  of  the 
number  of  panels  (DEC  2060).  For  VAX  11-750,  CPU  times  increase  by  approximately  2.5 
times. 
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Figure  13:  Mode  2  displacement  contours  for  the  DDH  propeller  blade  for  finite  element 
and  surface  panel  added  mass  matrices. 
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Figure  14:  Mode  3  displacement  contours  for  the  DDH  propeller  blade  for  finite  element 
and  surface  panel  added  mass  matrices. 
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Figure  16:  Mode  1  displacement  contours  for  the  P4388  propeller  blade  for  finite  element 
and  surface  panel  added  mass  matrices. 
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Figure  17:  Mode  2  displacement  contours  for  the  P4388  propeller  blade  for  finite  element 
and  surface  panel  added  mass  matrices. 
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Figure  19:  Mode  4  displacement  contours  for  the  P4388  propeller  blade  for  finite  element 
and  surface  panel  added  mass  matrices. 


Appendix  A:  Simplifying  assumptions  and  boundary  condi¬ 
tions 


This  appendix  presents  in  more  detail  the  assumptions  inherent  in  the  derivations  of 
the  linearized  pressure  and  surface  velocity  terms,  and  the  free  surface  boundary  condition 
used  in  the  surface  panel  formulation. 


Pressure  relation 


(31) 


The  Navier-Stokes  equations  describe  the  fundamental  equations  of  motion  of  a  fluid.  If 
we  consider  the  velocity  of  a  body  through  a  stationary  fluid  to  be  equivalent  to  the  motion 
of  the  fluid  about  a  stationary  body,  the  Navier-Stokes  equation  in  the  x  direction  for  an 
ideal  fluid  is 

du*  ,  du*  „  du*  ,  du*  1  dp 
dt  dx  dy  dz  p  dx 

where  u*  =  U  +  u,  etc.,  and  U,V,W  represent  body  velocities  and  u,v,w  represent  fluid 
disturbance  velocities  in  the  vicinity  of  the  body.  For  an  irrotational  flow  field,  potential 
flow  theory  can  then  be  used  to  define  a  velocity  potential  function  $,  where 

9  =  4*oo  +  4>  (32) 

with 

4>oo  =  Ux  +  Vy  +  Wz  (33) 

Here,  4>  is  the  disturbance  potential  in  the  fluid,  and  the  disturbance  velocities  are  defined 
as  potential  gradients 

d4>  _  d<f>  _  d4> 
dx  dy  dz 

Substituting  the  velocity  definitions  into  equation  (31)  and  rearranging,  we  obtain 

|£  +  (U  +  u)|^  +  (V  +  v)^  +  (W  +  «;)^  =  -i|| 

Integrating  equation  (35)  with  respect  to  x,  we  obtain  a  pressure  relation  in  the  form 


u  = 


(34) 


(35) 


?*  +  u^  + 
at  ax 


vii. + waA + i  r(«i)i + (»♦)« + ,?*f 

dy  dz  2  [  3x  Kdy ’  + Kdz ’ 


P 

P 


(36) 


A  linearized  form  of  this  relation  for  a  body  moving  in  the  fluid  can  be  obtained  by 
neglecting  the  higher  order  terms.  In  addition,  for  a  stationary  body  (equivalent  to  zero  far 
field  velocities  U,V,W),  this  equation  reduces  to 


P  =  ~P 


d4> 

dt 


(37) 


which  is  the  form  used  in  the  added  mass  matrix  derivation  given  in  this  report. 


Body  surface  boundary  condition 

The  time  derivatives  involved  in  the  body  surface  velocity  boundary  condition  will  also 
contain  nonlinear  terms  when  rigid  body  motion  is  included.  The  displacements  of  the  body 
are  defined  by  6(x,y,z,t),  for  which  the  total  derivative  is 


d6  _  86  86  dx  86  dy  86  dz 
dt  8t  +  8x  dt  +  dy  dt  dz  dt 


(38) 


Substituting  the  velocity  components  for  the  component  time  derivatives  in  (38),  the  full 
body  surface  velocity  becomes 


d6  86  ,rr  .86  .86  .86 

«  -  j-, + + +  + ">*; 


(39) 


The  terms  involving  the  disturbance  velocities  are  commonly  neglected  on  the  basis  of  order 
of  magnitude.  Introducing  the  transformation  matrix  to  relate  global  to  normal  directions, 
and  recalling  equation  (3),  the  complete  expression  for  the  control  point  normal  velocities 

fn  =  !T]{ :« >  +  U  |T,]  §f  +  V  py  §f  +  IV  Pil  If  («) 

Those  terms  involving  rigid  body  velocities  U,  V,  W  in  the  pressure  and  surface  boundary 
conditions  (38)  and  (40)  lead  to  additional  terms  in  the  nodal  force  vector.  For  example, 
the  inclusion  of  a  forward  speed  component  U  in  the  case  of  a  ship  hull  has  been  considered 
in  Reference  3.  The  complete  expression  for  the  force  vector  in  that  case  becomes 

{Ff}  =  -p[T]T[A]{H}[L]-1{T\{6} 

-  PU  ([T]r  [A]  [V,]  \L]~l  [T]  +  [If  [A]  [If]  [L]~l  [T\  {6} 

+  pU*  ([Tf  [A]  [V*]  [Lp1  [T]  {6} 

where  [V*]  is  a  coefficient  matrix  relating  source  strength  to  control  point  velocity  in  the  i 
direction. 

As  can  be  seen  from  equation  (41),  the  inclusion  of  the  convective  terms  does  not  change 
the  form  of  the  added  mass  matrix.  The  second  and  third  terms  of  the  right-hand  side  can 
be  considered  added  damping  and  added  stiffness  terms  respectively.  Their  effect  on  the 
dynamic  response  of  the  structure  is  not  considered  in  Reference  3,  nor  in  the  current  report 
since  the  finite  element  system  used  in  the  analyses  undertaken  here  does  not  include  a 
capability  to  include  a  damping  matrix  in  the  above  form.  The  added  stiffness  effects  could 
be  included;  however,  these  terms  are  of  higher  order  in  6,  and  are  probably  negligible  at 
realistic  forward  speeds. 


(41) 


•J 
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Free  surface  boundary  condition 

The  definition  of  the  free  surface  boundary  condition  represents  perhaps  the  most  serious 
approximation  in  both  the  finite  element  and  surface  panel  formulation  of  the  fiuid/structure 
interaction  problem.  As  discussed  in  Section  2,  the  surface  panel  method  uses  imaging  of 
sources  of  opposite  or  similar  sign  across  the  plane  of  the  free  surface  to  define  a  zero 
potential  (pressure  relief)  free  surface,  or  a  zero  normal  velocity  free  surface  respectively. 
In  the  finite  element  method,  the  free  surface  nodes  are  normally  specified  as  points  of  zero 
potential,  providing  again  the  pressure  relief  surface.  These  both  represent  approximations 
to  the  true  free  surface  boundary  condition,  a  discussion  of  which  is  presented  below. 

For  a  free  surface  which  includes  gravity  waves,  the  equation  of  motion  of  the  free  surface 


is 


(42) 


dz  _  dr) 
dt  dt 

where  z  =  h  +  rj  and  17  is  the  wave  amplitude,  and  h  the  depth  of  fluid.  For  potential  flow, 
we  have  then  the  kinematic  surface  condition 


<h 

dt 


dz 


(43) 


The  pressure  relationship  from  the  Navier-Stokes  equation,  including  the  surface  wave  pres¬ 
sure,  is 

d<b  ,  „ 

P  =  -9~£  +  P9V  (44) 

and  at  the  free  surface,  where  p  =  0,  equations  (42),  (43)  and  (44)  can  be  combined  to  give 


<P<I>  d<fr  . 


(45) 


Expanding  the  total  derivative,  and  using  again  the  potential  flow  relations  (32)  and  (33) 
with 


U  =  Uoo  V=W  =  0 


we  obtain 


d2</> 

dt* 


+  2  U. 


d2+  1  rr>  f  0 

“dxdt  +  H^dx*  +9dz~° 


(46) 


(47) 


as  the  complete  boundary  condition.  This  condition  is  usually  linearized  in  the  sense  that 
it  is  applied  at  z  =  0,  rather  than  at  the  dynamic  free  surface  z  =  rj.  It  is  instructive  to 
non-dimensionalize  equation  (47),  using 


t* 

'  =  4>*  4>  teal  t  t  =  —  X—X*L  Z  =  Z*L 

U) 


(48) 
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Here,  L  is  a  characteristic  length  over  which  the  potential  varies,  and  u>  a  characteristic 
radial  frequency.  Substituting  these  relations  into  (47)  and  dropping  the  prime  notation, 
we  have 

u2Ld'<t>  ,  2 Uoow  9*6  .  Uld*<t>  ,  96  „ 

9^i  +  ^L9^+9^  =  0  (49) 

as  the  free  surface  boundary  condition  which  the  potential  function  6  must  satisfy  in  addi¬ 
tion  to  the  Laplace  equation  and  the  body  surface  boundary  condition.  A  simple  solution 
does  not  exist  for  this  problem,  but  under  certain  conditions,  specific  terms  in  equation  (49) 
will  be  dominant,  and  a  solution  may  be  found  which  satisfies  the  simplified  form  of  the 
boundary  condition.  Several  cases  are  considered  below. 

Consider  first  a  zero  frequency  limit,  in  which  case  (49)  reduces  to 

£7*  9*6  96  ,  , 

=  0  <5°> 

A  fundamental  solution  which  satisfies  this  condition  and  the  Laplace  equation  can  be 
written  in  the  form19,20 

° = i + ik  fy  r /(*'  ■ '■ 1 fe' ' 1 " 

in  which  k  is  the  wave  number  and  /  is  an  exponential  function.  The  extra  term  represents 
an  integration  of  the  wave-spectrum  over  the  complete  free  surface.  In  the  limit  of  zero 
forward  speed,  equation  (51)  can  be  shown  to  reduce  to 


G  =  -  +  ^  (52) 

r  r* 

which  is  exactly  the  condition  obtained  from  the  positive  imaging  system.  This  method 
thus  satisfies  exactly  the  boundary  condition 

g  =  0  (53) 

in  the  zero  frequency  limit.  As  £/<x>  oo,  equation  (51)  can  be  shown  to  reduce  to 

G=  (54) 

r  r* 

which  is  the  condition  obtained  from  the  negative  imaging  system.  The  boundary  condition 


i  •'!  t’l  i1!  l 


m 

& 


corresponding  to  a  pressure  relief  surface  is  thus  exactly  satisfied  in  this  case.  Summarizing, 
for  the  ui  =  0  limit,  the  low  Froude  number  problems  can  be  approximated  with  the  pos¬ 
itive  imaging  method,  and  the  higher  Froude  number  problems  with  the  negative  imaging 
method. 

In  the  case  where  =  0  and  u  ^  0,  the  free  surface  boundary  condition  becomes 

fiMH  <-> 

Following  the  lines  of  the  previous  discussion,  it  can  be  seen  that  as  w  — ►  0,  the  boundary 
condition  can  be  satisfied  exactly  by  the  positive  imaging  method.  As  w  — ►  oo,  the  negative 
imaging  system  provides  a  correct  representation. 

The  ADMASS  system  currently  incorporates  both  the  positive  and  negative  imaging 
approach  for  the  free  surface  approximation,  and  the  above  discussion  indicates  the  types 
of  problems  for  which  the  formulation  will  provide  reasonable  accuracy.  In  the  general  case 
in  which  neither  limit  is  appropriate,  the  accuracy  will  depend  on  the  relative  magnitudes 
of  the  r on-dimensional  terms  in  the  complete  form  of  equation  (49). 

As  an  example,  consider  the  magnitudes  of  the  various  terms  which  are  relevant  to 
an  added  mass  calculation  for  a  ship  moving  at  a  realistic  forward  speed.  Take  typical 
parameters  as  U  =  10  m/aec,  L  =  5  m,  (draft).  As  a  function  of  w,  equation  (49)  becomes 

O.51w»|f  +  2.O^|t  +  2.O40  +  g=O  (57) 

For  u  =  1.0  (/  =  0.16Hz),  all  the  terms  are  of  the  same  order  of  magnitude,  and  poor 
accuracy  could  be  expected  from  the  imaging  method.  An  order  of  magnitude  dominance  in 
the  frequency  dependent  term  relative  to  the  final  term  is  achieved  for  w  =  4.4  (/  =  0.7 Hz), 
and  this  would  represent  the  approximate  lower  frequency  limit  for  which  the  negative 
imaging  system  might  be  expected  to  provide  adequate  accuracy  for  these  parameters.  Since 
this  limit  is  below  the  hull  bending  vibration  frequencies  of  almost  all  ships,  the  negative 
imaging  method  will  provide  added  fluid  mass  estimates  valid  for  natural  frequency  or 
dynamic  response  analyses.  The  positive  imaging  approach  appears  to  have  little  practical 
application  to  the  computation  of  added  fluid  mass  for  vibrating  structures  as  a  free  surface 
representation;  the  negative  system  has  provided  better  agreement  with  other  methods  in 
all  cases  considered  to  date. 

The  incorporation  of  the  more  complex  fundamental  solution  in  the  evaluation  of  the 
influence  coefficients  in  the  ADMASS  system  would  involve  a  combination  of  both  the 
current  imaging  system  and  an  evaluation  of  the  integral  term  of  equation  (51).  This 
approach  will  necessitate  the  introduction  of  a  numerical  integration  algorithm,  and  will 
add  to  the  computational  demand  in  that  section  of  the  system.  The  increased  demand 
may  be  justified  in  cases  where  the  limiting  conditions  satisfied  by  the  imaging  methods 
alone  are  not  appropriate. 
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Appendix  B  -  Surface  Panel  Integral  Evaluation 

In  this  appendix  it  is  shown  that  the  potential  of  a  planar  panel  with  uniform  source 
distribution  can  be  written  in  terms  of  expressions  for  the  velocities  induced  by  the  same 
panel.  The  velocity  relations  were  originally  derived  by  Hess  and  Smith.2 

The  number  of  edges  of  the  panel  will  be  denoted  Ne.  Let  xk,  k  =  1 ,Ne  denote  the 
the  corner  points  of  the  planar  panel  in  clockwise  order.  A  coordinate  system  is  chosen  so 
that  the  panel  lies  in  the  zy-plane;  hence,  z*  =  0. 

The  potential  at  x  =  (x,  y,  z)  induced  by  a  uniform  unit  source  on  the  panel  is 

<58> 

where  Tj  denotes  the  panel  surface.  The  associated  velocity  integral, 

VW=1,RF^  <59> 

has  been  solved  exactly  by  Hess  and  Smith2 

N. 

VZ  =  Y1V^  (60) 

*=i 

Vv  =  T,Vvk  (61) 

k-i 


where 


/ "»*,fc+icfc  “ 

•rctaa  "  h>+') 

V  zr*  / 

\  zrk+i  ) 

Vxt  =  (yt+1  -  y*)ln(H*)/dt,t+i 
Vvk  =  ~{xk+i  -  Xk)\n(Rk)/dk,k+i 
Rk  —  r*  +  r*+ 1  ~  dk,k+l 
*k  +  *■*+!  +  dk,k+l 

*k  =  y/(*~  xk)2  +  (y  -  yt)2  +  z2 
e*  =  (i  -  xk)2  +  z2 
hk  =  (y-y  k){x-xk) 

.....  Vfc+i  -  y* 

mk,k+i  - 

Zfc+i  -  Xk 

dk,k+ i  =  J{xk+\  -  xk)2  +  (y*+i  -  y*)2 
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In  these  formulae  it  is  understood  that  a  subscript  Ne  +  1  on  one  of  the  corner  points  is 
equivalent  to  a  subscript  1. 

Consider  the  integral  for  the  z-component  of  the  velocity: 


Let  V  be  the  two  dimensional  gradient  operator  with  respect  to  the  primed  coordinates. 

V  =  x—  +  y—  (75 

dxf^Vdy>  1 


V'(X)  =  /r, V  •  [f^|]  dxVs'’ 


which,  using  the  divergence  theorem,  can  be  reduced  to  a  line  integral  counterclockwise 
around  the  perimeter  of  the  panel. 


Vz  =  f  T^-^ds 
Jo r.-  x  -x' 


where  dTj  denotes  the  perimeter  of  the  panel,  ds  is  an  increment  of  arclength,  and  n  is  an 
outward  pointing  unit  normal  to  the  perimeter. 

The  edge  of  the  panel  between  x*  and  x*+i  can  be  parameterized  by 

x'  =  x*  +  (x*  -  x*+1)t  fort  €[-1,0]  (75) 

Note  that  the  vertices  must  be  passed  in  reverse  order  since  the  integral  proceeds  coun¬ 
terclockwise  around  the  perimeter  of  the  panel.  The  increment  of  arclength  is  then  da  = 
|xfc+i  -  Xjk| dt  =  dk,k+idt.  The  normal  is  given  by 

...  ...  (**  -  Xk+i)  x  * 

n  “ - 1 -  l76) 

Therefore,  n  •  x  =  ( y*  -  y*+i)dt.  Thus, 


f°  1 

=  -E(vk+i-y*)/_i|rr^* 


Comparison  with  equations  (60)  and  (63)  shows  that 


Now  consider  the  potential  of  equation  (58).  For  reasons  that  will  become  clear  shortly, 
rewrite  the  integral  as  follows. 


, ,  n  /*  z2  +  |x  -  x'|2  ,  ,  ,  ,  f  z2  ,  , 

*<x)  “  4  Tx^T~  y  ~  k  dy 

z2  +  |x-x'|2  ,  , 


ir,  |x-x'|s 


-dx'dy'  -  zVt 


(79) 


Now 


z2  +  |x  -  xf|2  _  2z2  +  (s  -  X *)2  +  (y  -  y')2 

|x  -  x'is  “  |x  x'l3 

l(z-x')£+(y-y')y2 


=  -V 


x  -  x' 


(80) 


Using  the  divergence  theorem  again, 


*(x)  =  -  L, 


n  •  [(x  -  x')x  +  (y  -  y')y] 


Ix-xM 


ds  -  zVt 


_  _  f 

'  Jar,  |x  -  x'| 


(81) 


(82) 


But,  parameterizing  the  sides  of  the  panel  as  before,  one  obtains 
n  •  (x  -  x')ds  =  [(xjt  -  x*+i)  X  z)  •  (x  -  x')dt 
=  I(x  -  x')  x  (Xfc  -  xfc+i)|  •  zdt 
=  [(x  -  x*  -  (xjfc  -  xt+i)t)  x  (xt  -  x*+i)]  •  zdt 
=  [(x  -  xt)  x  (x*  -  Xfc+i)]  •  zdt 
=  [(*  -  **)(y*  -  y*+i)  -  (y  -  y k){xk  -  xk+i))dt 
Therefore,  making  use  of  equations  (78),  (81),  and  (82), 

y0  | 

<£(x)  =  +  £[(*- z*)(yfc-u  -  yiO  -  (y  -  y*)(zM^  -  ^)]  / 

k=i  |X  X| 

=  ~zV*  -  ~  z*)(y*+i  -  y*)  -  (y  -  y*)(x*+ 1  -  xk)]^Rk^ 

t=i  dk-k+1 

N. 

=  -zVa  -  22(x-  xk)Vzk  +  (y  -  yk)Vvk  (83) 


:dt 


k=l 


Hence,  the  potential  integral  can  be  written  in  terms  of  the  same  algebraic  expressions  as 
the  velocity  integral.  The  extra  computing  time  needed  to  calculate  the  potential  integral 
in  addition  to  the  velocity  integrals  is  very  small. 
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Multipole  Expansion  for  Large  Separation 

When  the  separation,  r,  «f  the  panel  and  the  field  P»i«  ~neid„^ie  !»*>.< «.»- 
tione  (58)  and  (59)  can  be  expanded  in  terms  of  the  email  quan  y  |  I/I 
expression  for  the  potential  may  be  written  as 


E  4-»s(;)+  E  E  + 

r  p€(^»)  dpKJ 


and  for  the  velocity  components  as 


E  E.|«’5=b(;)+-1 


«€(*,!/)  *€(*.!/) 


for  p  6  (x,y,«) 


where 


/(°)  =  f  dx'dy' 

Jr, 

iff  =  [  x' dx'dy'  ;  4X)  =  /  y' dx'dy1 

JVj  Jr> 

iff  -  f  x'2  dx'dy'  ;  iff  =  j^x'y' dx'dy'  ;  4?  = 

JTj  i 


and  so  on.  The  derivatives  of  l/r  are 


pe(w) 

&(;)~*+£  »,«<«..> 

3s  /1\  3p5g,  ,  3ggg,  ^  SsSpg  _  15pgs  p  9)Se  (x,y,z) 

+  r&  '7 


e  .  .1  delta  Note  that  the  panel  moments  are  properties  of  the 

:S^n“  t, “ie“d  n^n  the  point  at  which  the  .elocit,  i.  to  he  evaluated 
Details  of  the  method  hy  which  they  are  evaluated  m  the  ADMASS  V***™*”’ 


Details  of  the  methoa  Dy  wmcn  - - - ;  '  ...  thp  Danei 

Hally.*1  The  origin  of  the  panel  coordinates  has  been  chosen  to  be  the  centroid  of  the  P 

so  that  iff  =  iff  =  0;  hence,  the  second  term  in  the  expansion  vanishes.  This  choice  also 


IVaVhV.V.'j 


ensures  that  the  multipole  expansion  can  be  used  whenever  r  is  much  larger  than  the  mean 
dimensions  of  the  panel.  Up  to  third  order, 


and 


/(°)  (3x2  -  r')liV  ZzyllV  (3 y2  -  r2)/$ 

+ - 2r® - +  + - 2r® - 


Tr  x/<°>  3x(3r2  —  5xJ)/*x^  3y(r2  -  5x2)J$  3x(r2  -  5y2)I$ 

y  &  — - - - - 


2r7 


2r7 


Vy  w 


y 3y(r2  -  5x2)/|^  3x(r2  -  5y2)/ffi  3y(3r2  -  5y2)/^ 


2r7 


2r7 


z/(°)  3z(r2  -  5x2)/i*^  15xyz/i^  3z(r2  -  5y2)2^ 

*  "  r3  2r7  r7  2 ? 


(92) 

(93) 

(94) 

(95) 


In  the  ADMASS  system,  the  exact  solutions  corresponding  to  equations  (60),  (61),  (62) 
and  (83)  are  used  when  r  is  small.  At  moderate  separations,  the  first  three  term  expansions 
of  the  expansions  of  equations  (93)-(95)  are  utilized.  At  large  separations,  the  first  term  of 
those  expansions  is  utilized  in  the  integral  evaluation. 


Appendix  C  -  Program  Execution 

This  appendix  describes  the  input  requirements  and  program  prompts  in  the  ADMASS 
system.  The  system  can  be  run  interactively  as  a  series  of  programs  ADMS1  -  ADMS4, 
or  as  a  single  system  from  a  command  file  in  batch  mode.  The  interactive  form  will  be 
discussed  here:  the  necessary  responses  are  the  same  in  either  approach. 

The  ADMASS  system  begins  with  the  ADMS1  routine,  which  accesses  the  structural 
and  fluid  geometry  files  PREFIX.GOM  and  PREFIX.  AMD.  The  format  of  these  files  should 
be  exactly  as  described  in  Reference  5,  with  the  exception  that  the  first  element  group  in 
the  .AMD  file  must  contain  all  of  the  interface  elements.  Mixing  of  interface  element  types 
is  not  permitted  at  the  present  time,  and  only  element  types  1  and  3  are  permitted.  The 
.GOM  file  need  only  contain  the  structural  nodal  coordinates  since  element  definitions  are 
not  read.  The  ADMS1  routine  rarely  takes  more  than  a  few  seconds  of  CPU  time  for  com¬ 
pletion.  The  prompts  and  expected  responses  (in  brackets  -  use  capitals)  for  this  routine 
are  as  follows. 

Enter  5  character  prefix  of  VAST  files  (Character) 

Reverse  the  interface  node  order?  (Y,N) 

If  the  response  to  the  second  prompt  is  Y,  the  corner  node  numbering  of  all  panels  will 
be  reversed  from  that  defined  in  the  .AMD  file.  This  causes  a  reversal  in  the  panel  normal 
orientation  in  the  event  that  connectivity  has  been  input  consistently  in  reverse  order.  To 
avoid  the  necessity  of  this  option,  panels  should  be  defined  in  a  counter-clockwise  sense 
viewed  from  the  fluid. 

Are  interface  elements  uniquely  defined?  (Y,N) 

This  option  is  necessary  when  the  thin  body  assumption  has  been  used  in  the  interface 
element  definitions  and  a  second  surface  needs  to  be  generated.  If  all  panels  would  corre¬ 
spond  to  a  unique  fluid  element,  this  option  is  not  required.  If  the  panels  are  not  unique, 
several  options  will  be  presented: 

Are  solid  elements  used  in  structure?  (Y,N) 

If  20-node  solid  elements  have  been  used  to  model  the  structure,  the  second  surface  can 
be  defined  by  specifying  node  number  increments  for  the  corner  and  midside  nodes  of  the 
available  surface.  These  increments  will  require  some  knowledge  of  the  model  generation 
system,  and  the  validity  of  the  second  surface  is  obviously  dependent  on  the  compatibility 


of  the  node  numbering  system  with  this  algorithm.  The  prompts  for  this  option  are 


Enter  corner  node  increment  for  surface  2  (Integer) 

Enter  midside  node  increment  for  surface  2  (Integer) 

One  other  option  is  available  for  generation  of  the  second  surface  in  the  case  that  shell 
element  modelling  option  1  has  been  used  in  the  structural  model  (see  Reference  1  pp.  02- 
El. 4).  This  option  appears  if  N  has  been  entered  in  response  to  the  solid  element  prompt 
above. 


Shell  element  option  1  used?  (Y,N) 

If  Y  is  entered  here,  the  program  automatically  uses  the  parameter  NDN  (number  of  dis¬ 
placement  nodes)  defined  in  the  .GOM  file  as  the  nodal  increment  for  all  surface  1  nodes. 
This  again  is  not  a  generally  applicable  approach,  but  has  been  developed  with  the  solid 
modelling  option  to  be  compatible  with  current  versions  of  the  propeller  modelling  system.15 
If  neither  solid  nor  shell  element  option  1  has  been  used  in  the  model,  a  second  surface  must 
be  defined  by  alternative  means. 

In  all  cases  where  the  8-node  interface  element  has  been  used,  two  refinement  options 
are  available: 

Edge  1  panel  refinement  required?  (Y,N) 

Edge  2  panel  refinement  required?  (Y,N) 

These  options  divide  each  panel  defined  in  the  .AMD  file  into  two  elements,  using  the 
midside  nodes  of  two  opposite  edges  and  the  4  corner  nodes.  Only  one  of  these  options  can 
be  used  at  a  time,  since  an  additional  node  within  the  interface  element  would  be  required 
for  refinement  in  both  directions.  Edge  1  of  an  element  is  defined  by  the  first  2  nodes  in  the 
interface  element  definition;  edge  2  is  defined  by  the  second  and  third  nodes  (see  Figure  5). 
The  program  should  now  return  the  following  information: 

Fluid  geometry  file  read  completed. 

Number  of  panels  is: _ 

Number  of  interface  nodes  is: _ 


These  values  should  be  checked  against  expectations  to  determine  if  the  refinement  al¬ 
gorithms  are  working  properly. 

A  header  for  identifying  output  disc  files  will  now  be  required, 


Enter  a  title  for  identifying  output  files  (Character) 

Since  all  file  names  are  fixed  within  the  ADMASS  system  and  the  files  are  in  binary  form, 
the  header  serves  as  the  only  check  in  subsequent  routines  that  the  correct  files  are  being 
used. 

Panel  properties  are  now  calculated,  and  the  following  lines  should  appear  as  calcula¬ 
tions  proceed 

Calculating  panel  properties... 

Panel  Property  and  Transform  files  written. 

Exiting  ADMSl.. 

The  output  files  are  titled  PANPRP.MAS  and  TMTRX.MAS  respectively.  Several  error 
checks  are  made  during  the  panel  property  calculations.  The  most  common  error  is  the 
nodal  connectivity  definition,  which  can  usually  be  corrected  with  the  node  reversal  option. 
If  errors  in  element  connectivity  result  in  non-quadrilateral  shapes,  the  error  message 

The  control  point  of  panel  (  )  does  not  lie  within  the  panel  boundaries 

will  usually  be  triggered.  To  facilitate  debugging  of  interface  grids,  the  ADMSl  routine 
writes  the  interface  element  nodes  (by  element),  the  panels  a  particular  node  is  connected 
to  (by  node)  and  the  panel  corner  coordinates,  area,  and  normal  components  to  a  formatted 
disc  file  ADMS1.DAT. 

The  ADMS2  program  can  be  run  on  successful  completion  of  the  ADMSl  routine.  This 
program  calculates  the  influence  coefficients  relating  the  panel  source  strength  to  control 
point  (panel  centroid)  velocity  potential  and  normal  velocity.  This  system  is  quite  efficient, 
but  can  take  several  minutes  of  CPU  time  for  a  larger  problem  (>  200  panels).  The  follow¬ 
ing  information  and  prompts  should  appear: 

Entering  ADMS2... 

Reading  panel  properties  from  PANPRP.MAS  ...  Done 
Header  is: 

Prefix - Time,  Date 

Enter  free  surface  boundary  condition 
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1.  1/R  +  1/R’  :  Symmetry  reflection  of  structure  across  surface. 

2.  1/R  -  1/R’  :  Pressure  release  boundary. 

3.  No  free  surface. 

The  response  to  this  prompt  is  dependent  on  the  particular  problem.  In  all  cases  con¬ 
sidered  to  date,  option  2  has  provided  a  better  approximation  to  the  free  surface  than 
option  1.  If  the  body  is  completely  submerged,  option  3  is  required.  The  next  prompt, 

Calculate  influence  coefficients  for  velocity  components?  (Y,N) 

defines  influence  coefficient  matrices  unrelated  to  the  added  mass  problem,  and  a  N  re¬ 
sponse  should  be  entered  here.  The  symmetry  option  may  be  used, 

X-Z  plane  of  symmetry?  (Y,N) 

for  half  models,  noting  again  the  specific  coordinate  system  that  is  used  within  the  program. 
A  further  prompt  concerning  antisymmetric  influence  coefficient  matrices 

Calculate  antisymmetric  influence  coefficients?  (Y,N) 

can  also  be  answered  N  for  an  added  mass  calculation.  The  ADMS2  program  then  begins 
the  calculation  of  the  influence  coefficient  matrices,  informing  the  user  of  progress  on  an 
ongoing  basis.  The  binary  output  files  of  this  routine  are  INFPHI.MAS  and  INFCFF.MAS. 

The  ADMS3  program  inverts  the  matrix  stored  on  the  disc  file  INFCFF.MAS,  producing 
the  file  IN  VERS. MAS.  This  is  generally  the  lengthiest  section  of  the  added  mass  matrix 
generation.  The  minimum  CPU  time  demand  of  the  current  algorithm  is  proportional  to 
between  (NP)2  and  (NP)3,  where  NP  is  the  number  of  surface  panels.  The  actual  CPU  time 
will  depend  on  the  convergence  performance  for  the  matrix  being  inverted.  For  models  with 
fluid  on  one  side  only,  the  CPU  time  is  usually  quite  close  to  the  minimum. 

The  method  uses  Gauss  elimination  and  back  substitution  in  the  solution  of  blocks  along 
the  matrix  diagonal;  the  number  of  blocks  is  the  only  input  to  this  routine.  In  test  runs, 
an  optimum  efficiency  appears  to  be  obtained  when  the  block  Bize  is  10  to  15  percent  of  the 
matrix  dimension.  The  solutions  of  these  blocks  are  then  used  in  a  Seidel  form  of  iterative 
solution,  using  the  identity  vectors  as  right-hand  sides. 

The  ADMS3  routine  will  print  an  iteration  count  when  the  number  of  iterations  for 
convergence  of  any  vector  exceeds  five.  If  poor  convergence  performance  is  indicated  by 
the  routine,  a  restart  with  a  smaller  number  of  blocks  is  recommended.  The  maximum 
number  of  iterations  performed  within  the  routine  is  twenty;  if  this  count  is  indicated  in 
the  statistics,  the  run  should  be  aborted. 


The  following  format  should  appear  when  this  program  is  run: 

Entering  ADMS3.... 

Total  number  of  panels  = _ 

Enter  number  of  blocks  desired  for  decomposition 

(Maximum  is  40  -  0  for  default)  (Integer) 

The  number  of  blocks  need  not  be  a  round  divisor  of  the  matrix  size  NP.  On  entry  of 
the  desired  number  of  blocks,  the  program  informs  the  user  of  progress  on  an  ongoing  basis. 

If  the  iteration  count  for  any  vector  exceeds  five,  the  following  message  will  be  printed: 

Poor  convergence  performance  for  this  matrix:  _ iterations  required 

The  tolerance  sought  in  the  iteration  is  less  than  0.001  in  the  current  version.  This  error 
limit  has  little  effect  on  the  accuracy  of  the  added  mass  matrix  for  well  conditioned  coef¬ 
ficient  matrices,  but  can  drastically  increase  the  CPU  time  required  for  matrix  inversion. 
The  current  setting  is  considered  to  be  a  good  balance  between  accuracy  and  computational 
efficiency,  and  thus  has  not  been  made  an  interactive  variable. 

The  final  routine  in  the  ADMASS  calculation  suite  is  the  ADMS4  program,  which 
performs  a  series  of  matrix  multiplications  to  produce  the  fluid  added  mass  matrix.  As 
discussed  in  Section  2  of  this  report,  the  source  and  dipole  formulations  of  the  surface  panel 
method  represent  different  multiplication  orders,  and  the  option  for  one  or  the  other  is  given 
here.  This  routine  can  be  reasonably  CPU  intensive,  often  taking  about  2/3  the  CPU  time 
used  by  the  ADMS3  program.  The  inputs  to  this  routine  are  the  panel  formulation  desired 
and  the  density  of  the  fluid.  The  program  uses  the  files  TMTRX.MAS,  INFPHI.MAS,  and 
INVERS.MAS;  again,  progress  statements  are  given  as  operations  are  begun  and  completed. 
The  output  of  this  routine  is  the  fluid  mass  matrix  PREFIX.T36,  compatible  with  the  VAST 
finite  element  analysis  system. 

Post-processing  options  to  read  selected  entries  in  the  added  mass  matrix,  compute 
rigid  body  masses,  condense  or  to  symmetrize  the  matrix  have  been  incorporated  in  a  pro¬ 
gram  ADMSP.  For  matrix  condensation,  the  program  accesses  the  TMTRX.MAS  file  for 
information  on  the  number  of  displacement  nodes.  A  warning  is  given  to  ensure  that  the 
TMTRX.MAS  matrix  corresponds  to  the  added  mass  matrix  being  reduced.  The  prompts 
for  this  program  are 

Enter  name  of  T36  file  to  read  (5  character  prefix) 


Enter  post-processing  option  required 


1:  Read  selected  entries  in  the  matrix 
2:  Rigid  body  mass  additions 
3:  Reduction  to  displacement  nodes 
4:  Symmetrize  mass  matrix 
5:  Stop 

In  cases  where  both  a  reduction  and  symmetrization  are  desired,  the  reduction  should  be 
done  first,  and  this  menu  will  be  offered  on  completion  of  the  reduction.  The  program  will 
allow  only  one  reduction  of  the  mass  matrix.  The  reduction  process  is  quite  CPU  intensive, 
and  will  take  up  to  about  20  percent  of  the  CPU  time  required  to  create  the  mass  matrix. 
The  symmetry  option  is  not  CPU  demanding. 


References 


1.  ‘Studies  on  Fluid  Structural  Vibration  Using  Boundary  Element  Method’,  Mitsubishi 
Technical  Bulletin  No.  172,  June,  1986. 

2.  Hess  A.M.O.,  Smith  J.L.,  ‘Calcualtion  of  Non-lifting  Potential  Flow  About  Arbitrary 
Three-dimensional  Bodies’,  Report  E.S.  40622,  Douglas  Aircraft  Co.,  1962. 

3.  Vorus  W.S.,  Hylarides  S.,  ‘Hydrodynamic  Added-Mass  Matrix  of  Vibrating  Ship 
Based  on  a  Distribution  of  Hull  Surface  Sources’,  Trans.  SNAME,  Vol.  89,  1981, 
pp.  397-416. 

4.  Hylarides  S.,  Vorus,  W.S.,  ‘The  Added  Mass  Matrix  in  Ship  Vibration  Using  a  Source 
Distribution  Related  to  the  Finite  Element  Grid  of  the  Ship  Structure’,  Int.  Ship¬ 
building  Progress,  Vol.  29,  No.  330,  1982. 

5.  ‘Vibration  and  Strength  Analysis  Program  VAST,  Version  4,  User  Manual’,  DREA 
Contract  Report  CR86/429,  Martec  Ltd.  1986. 

6.  Lee  J.,  ‘A  Potential  Based  Panel  Method  for  the  Analysis  of  Marine  Propeller  Blades 
in  Steady  Flow’,  Report  87-13,  Massachusetts  Institute  of  Technology,  Department  of 
Ocean  Engineering,  Cambridge,  Mass.,  1987. 

7.  Norwood  M.,  'Application  of  the  Finite  Element  Method  for  the  Calculation  of  the 
Hydrodynamic  Added  Mass  of  Marine  Propeller  Blades’,  DREA  Informal  Communi¬ 
cation,  1977. 

8.  Hally  D.,  ‘POTFLO  -  A  Suite  of  Programs  for  Calculating  Potential  Flow  About  Ship 
Hulls’,  DREA  Technical  Memorandum,  in  review. 

9.  Sarpkaya,  T.,  Isaacson,  M.,  Mechanics  of  Wave  Forces  on  Offshore  Structures.  Van 
Nostrand  Reinhold,  New  York,  1981. 

10.  Lewis,  F.  M.,  ‘The  Inertia  of  Water  Surrounding  a  Vibrating  Ship’,  Trans.  SNAME, 
Vol.  37,  1929,  pp.  1-20. 

11.  Hicks,  A.  N.,  ‘The  Theory  of  Explosion  Induced  Ship  Whipping  Motions’,  Naval 
Construction  Research  Establishment  Report  R579,  Dunfermline,  U.K.,  1972. 

12.  Hally  D.,  ‘Numerical  Approximation  of  Ship  Hulls’,  DREA  Technical  Memorandum 
TM85/207, 1985. 

13.  Hally  D.,  ‘Representation  of  Hulls  with  Knuckles’,  DREA  Technical  Memorandum 
TM86/209, 1986. 


HBWBBIWff 


14.  Glenwright  D.,  Hutton  S.,  ‘Experimental  and  Finite  Element  Investigation  of  Added 
Mass  Effects  on  Ship  Structures’,  DREA  Contract  Report  CR87/464, 1987. 

15.  ‘Vibration  and  Stress  Analysis  of  Marine  Propellers  by  Finite  Element  Methods’, 
DREA  Contract  Report,  Martec  Ltd.,  1983. 

16.  Smith  D.R.,  Slater  J.E.,  DREA  RN/H/79/4,  Informal  Communication. 

17.  Slater  J.E.,  Smith  D.R.,  ‘Prediction  of  Static  Stress,  Natural  Frequencies  and  Mode 
Shapes  for  the  DTNSRDC  4388  Highly  Skewed  Propeller  Blade  by  Finite  Element 
Analysis’,  DREA  RN/H/82/19,  Informal  Communication. 

18.  Dhir  S.K.,  Sikora  J.P.,*  Holographic  Study  of  Vibration  Characteristics  of  a  Highly 
Skewed  Propeller  Blade  Model’,  DTNSRDC  Technical  Note  193,  August,  1971. 

19.  Doctors  J.L.,  Beck  R.F.,  ‘Numerical  Aspects  of  the  Neumann-Kelvin  Problem’,  J.  of 
Ship  Research,  Vol.  31,  No.  1,  March  1987,  pp.  1-13. 

20.  Wehausen  J.V.,  Laitone,  E.V.,  ‘Surface  Waves,’  Encyclopedia  of  Physics,  Vol.  9, 
Flugge,  Ed.  Springer  Verlag,  Berlin, 1960. 

21.  Hally  D.,  ‘Calculation  of  the  Moments  of  Polygons’,  DREA  Technical  Memorandum 
TM87/209,  1987. 


_ Unclassified 

StCOhfTY  CLASSIFICATION  OF  FOKM 
(highest  elessificetien  of  Title,  Abstract.  Keyword!! 


DOCUMENT  CONTROL  DATA 

(Security  clMsdication  of  till*,  body  of  abstract  and  indexing  annotation  must  ba  antorad  wHan  tha  overall  document  is  classified) 


1,  ORIGINATOR  tote  name  and  address  of  the  or  sanitation  preparing  the  document 
Organisations  tar  whom  the  document  was  prepared,  e.p  Establishment  sponsoring 
•  contractor's  report  or  tasking  agency,  ara  tnttrtd  in  section  LI 

Defence  Research  Establishment  Atlantic 


2.  SECURITY  CLASSIFICATION 

(overall  security  classification  of  the  document, 
including  special  warning  terms  li  applicable) 

Unclassified 


3.  TITLE  (the  complete  document  title  as  indicated  on  the  title  pege.  Its  clissificetion  should  be  mdiceted  by  the  appropriate 
abbreviation  {S.C.R  or  U)  in  parentheses  efter  the  title.) 

A  SURFACE  PANEL  METHOD  FOR  THE  CALCULATION  OF  ADDED  MASS  MATRICES  FOR 
FINITE  ELEMENT  MODELS 


4.  AUTHORS  (Last  name,  first  name,  middle  initial.  If  military,  show  rank.  e.g  Doe.  Met  John  E.) 


Vernon.  T.A..  Bara,  B.«  Hally,  D. 


5.  DATE  OF  PUBLICATION  (month  and  year  of  publication  of 

6a  NO.  OF  PAGES  hotel 

6b.  NO.  OF  REFS  (total  cited  in 

document) 

containing  information.  Include 

document) 

February  1988 

Annexes,  ^pendices,  etc.) 

21 

6.  DESCRIPTIVE  NOTES  (the  category  of  the  document,  e.g  technical  report,  technical  note  or  memorandum.  If  appropriate,  enter  the  type  of 
report,  e.g.  inte'iir,  progress,  summary,  annual  or  final.  Give  the  inclusive  dates  when  a  specific  reporting  period  >s  covered.) 


DREA  Technical  Memorandum _ 

8.  SPONSORING  ACTIVITY  (the  name  ef  the  deportment  protect  office  or  laberotory  sponsoring  the  research  and  development.  Include  the 
address.) 


9a  PROJECT  OR  GRANT  NO.  (i(  appropriate,  the  applicable  reseerch 
and  development  project  or  grant  number  under  which  the  document 
was  written.  Please  specify  whether  project  or  grant) 


9b.  CONTRACT  NO.  (if  appropriate,  the  applicable  number  under 
which  the  document  was  written) 


10a.  ORIGINATOR'S  DOCUMENT  NUMBER  (the  official  document 
number  by  which  the  document  is  identified  by  the  originating 
activity.  This  number  must  be  unique  to  this  document.) 

DREA  TECH.  MEMORANDUM  88/203 


10b.  OTHER  DOCUMENT  NOS.  (Any  other  numbers  which  may 
be  OSSipned  this  document  either  by  the  originator  or  by  the 
Spenser) 


1 1.  DOCUMENT  AVAILABILITY  (ary  limitations  on  further  dissemination  of  the  document,  other  then  these  imposed  by  security  classification) 


J  Unlimited  distribution 

I  Distribution  limited  (.  defence  deportments  aid  defence  contractors;  further  distribution  only  is  ipproved 

)  Distribution  limited  to  defence  departments  aid  Conodion  defence  csntrectors;  further  distribution  only  os  ^proved 

)  Distribution  limited  to  government  depertments  aid  agencies;  further  distribution  only  os  approved 

)  Distribution  limited  to  defence  depatments;  further  distribution  only  os  vproved 

)  Other  fpleoso  specify): 


12.  DOCUMENT  ANNOUNCEMENT  (ary  limitation  to  the  bibliogrmhic  atnouncement  of  this  document  This  will  normally  correspond  to 
the  Oacumom  Avoiiabifty  (111.  Hewovor,  whore  further  distribution  (beyond  the  audience  specified  in  11)  is  possible,  o  wider 
awounctmom  audience  may  be  selected.) 

Unlimited 


Unclassified 


StCUhlTV  CLASSIFICATION  OF  FOKM 


DCD03  2/06/87 


_ Unclassified 

SfCUAITV  CLASSIFICATION  Of  FOAM 


13  ABSTRACT  (  i  brief  and  factual  summary  of  the  document  It  may  also  appear  elsewhere  *r  body  of  the  document  itself.  It  is  highly 
da  sir  ip  la  that  the  ^street  of  classified  documents  be  unclassified.  Each  paragraph  of  the  ebsti  r  ■  thali  begin  with  an  indication  of  the 
security  classification  of  the  information  m  the  peregrin  (unless  the  document  itself  is  under  J  represented  as  IS).  1C).  Iff),  or  (U) 

It  is  net  necessary  ta  include  here  Attracts  in  both  office!  languages  unless  the  test  is  bilingum,. 


A  method  of  generating  the  added  fluid  mass  matrix  for  use  in  dynamic 
analyses  of  submerged  structures  is  presented.  This  method  uses  a 
distribution  of  panel  singularities  in  the  form  of  sources  or  dipoles  on 
the  fluid/ structure  interface  to  represent  the  velocity  potential  in  the 
fluid  surrounding  the  structure.  The  fluid  added  mass  matrix  is  calcula¬ 
ted  by  relating  the  pressure  field  in  the  fluid  to  the  structure  surface 
accelerations  via  the  Euler  equation  and  the  assumptions  of  potential  flow 
theory.  Based  on  the  results  of  a  variety  of  investigations  of  computat¬ 
ional  performance,  convergence  and  dynamic  response  involving  floating 
cylinders,  ship  hulls,  and  propeller  blades,  the  surface  panel  method  is 
shown  to  provide  an  effective  alternative  to  the  finite  element  added  mass 
matrix  system.  The  surface  panel  program  suite  has  been  developed  in  a 
form  compatible  with  an  existing  finite  element  analysis  package,  although 
the  dependence  on  the  particular  solution  system  is  not  large. 

,  'pot&hK*  i  fric,';JzS 

/  '  x  J'  ■  , ") 


Cl  -T. 


^  T  5  V  1 


::r.?5S 

l.  y- 

J  J  lie 


14.  KEYWORDS.  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a  document  end  could  be 
helpful  in  cataloguing  the  document  They  should  be  selected  so  Ihol  no  security  classification  is  required  Identifiers,  such  os  equipment 
model  designation,  trode  name,  military  protect  code  name,  geographic  location  may  also  be  included.  II  possible  keywords  should  be  selected 
ftom  o  published  thesaurus,  eg  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  is  not  possible  to 
select  indexing  terms  which  ore  Unclassified,  the  classification  of  each  should  bt  indicated  as  with  the  title) 


Mass  matrix 
Added  mass 

Fluid/structure  interaction 

Interaction 

Dynamics 

Finite  element 

Panel 

Surface  panel 
Boundary  integral 
Potential  flow 


Unclassified 


SfCUAITV  classification  of  foam 


